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POLYMERIC PRECISION 


This ebulliometer—being assembled the master glass- 
blower his own lathe our Central Research Laboratory. 

The ebulliometer uses two vapor lift pumps such those found 
electric coffee percolators. One circulates polymer solution 
while the other circulates boiling solvent—each liquid stream 
directed over temperature detector sensitive that difference 
one-millionth may sensed. This device allows the 
measurement high molecular weights related the physical 
properties polymers. 

Such C-I-L research contributes higher standards per- 
formance the plastics industry new and better plastic 
products and—through better technical service—to greater 
success the application these new materials modern living. 
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Treatment Uranium Leach Plant 
Produce High Purity Uranium 


Three liquid-liquid extraction processes were 
tested continuous equipment, bench 
determine their relative merits the purification 
uranium found ion exchange eluates various 
leach plants. Tri-n-butyl phosphate and 
dibutyl butylphosphonate were used extractants 
nitrate eluates. tertiary amine was used the 
case sulphuric acid eluates. Optimum conditions 
for high recovery and low reagent costs were estab- 
and the adaptability the processes exist- 
ing leach plants discussed. The suitability the 
ammonium diuranate product source ceramic- 
grade UO, was also investigated. Suggested changes 
the adopted flowsheets are given, for further 
pilot plant scale. 


Status Canadian Uranium Industry 
uranium industry now well established Canada and 
the processes used the production leach plant concen- 
and reactor-grade oxide metal are known the 
netallurgist. paper will refer times details the 
production flowsheets, and will useful this time define 
where the present work fits the uranium picture. 


All Canadian uranium leach plants, with the exception 
the Eldorado plant Beaverlodge, Saskatchewan, use sulphuric 
leaching followed ion exchange. (The Eldorado plant 
has adopted the sodium carbonate leach process). 
acid leach plant flowsheet showing the important 
steps the process shown Figure since the main topics 
this paper are associated with the last three operations, 
few comments these may interest. 

Recovery and concentration uranium from acid sulphate 
stationary moving-bed anion exchange resins 
both efficient process and cheap terms the overall cost. 
Resin life measured years, while basic chemicals used 
for only cents the cost per pound uranium 
when nitrate elution and moving-bed columns are 
Precipitation most plants two stages. the first stage, 
added 3.0 3.5 precipitates most the iron and 
and, the second stage, magnesia, caustic soda, 
ammonia recovers diuranate relatively free iron, thorium, 


‘Manuscript received March 27; accepted July 28; 1961. 

Metallurgy Division, Mines Branch, Department Mines and 
Technical Surveys, Ottawa, Ont. 

paper presented the Joint Chemical Engin- 
Conference, Cleveland, Ohio, May 7-10, 1961. 
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phosphates and halides. This second product dried and then 
shipped drums customs refinery exported other 
countries mainly the U.S.A. and Britain, and, minor 
quantities, Germany and Japan. 

dissolved nitric acid and the nitrate extracted with 
25% phosphate (TBP) kerosene; the extract 
scrubbed and finally stripped with water. This water strip, 
this point, either evaporated recover uranium trioxide 
neutralized precipitate ammonium diuranate (ADU). Depend- 
ing the demand for uranium hexafluoride feed the 
isotope separation plant for natural uranium metal, the trioxide 
processed accordingly. Recent the use 
sintered uranium dioxide reactor fuel have emphasized the 
ammonium diuranate route. 


alternate route for the production the hexafluoride 
from leach plant concentrate has been developed the U.S.A. 
whereby direct hydrofluorination the concentrate, followed 
fluorination and distillation, can produce suitable feed 
the isotope plant. 


Short-cuts Refined Uranium 
Looking back the leach plant, where operating costs are 
constantly being reduced improved maintenance, refine- 


OPERATION | REAGENTS EQUIPMENT PRODUCT 
| 1) Crush and Grind Steel \9/ Mine Ore 
| 2) Leach Sulphuric Acid 
Air 
Chlorate 
3) Filter | Flocculonts 
Clarify Lime 


| (for waste disposal) 


Exchange 


Nitric Acid 
or Salt— 
Sulphuric Acid 
Barren Effluent 
Adsorption (to waste) 


5) Precipitation jlime 


Ist Stage Lime and ammonia 

| Iron - Gypsum Cake 
(to leach) 

| 2nd Stage Ammonia 

| 

| 6) Drying and Fuel Uranium Cake 


Figure 1—Uranium acid leach plant flowsheet. 
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OPERATION | REAGENTS | FLOWSHEET 


PRODUCT 


Uranium 


1) Dissolve Nitric Acid 


2) Extract 25-30% TBP 


Scrub Nitric Acid 
Strip Water 


Nitric 
| Acid 
3) Evaporate | Fuel Recovery 
or 


Precipitate Ammonia 


Ammonium 


| 4) Reduce Cracked ammonia | Nitrate 
Press or | 


Sinter Hydrogen 


5) Hydrofiuorinate Hydrofluoric Acid | Ceramic 
ond UO, Pellets 


| Fluorinate Fluorine | UF 
| | to isotope 


| F. separation 
6) Reduce Magne: 2 

| ) gnesium | ———_—— U Metal 


Figure 2—Uranium refining and fabrication flowsheet. 


ment control, recycling some waste products, 
apparent that simplification can also achieved by-passing 
the precipitation, drying and packaging steps and refining the 
uranium directly from ion exchange eluates. The merging 
plant and refinery produce “reactor grade” ammonium 
diuranate was the object some laboratory and pilot plant 
studies the Mines Branch. The feasibility, cost and adapta- 
bility three liquid-liquid extraction processes, applied ion 
exchange eluates, were the main concern these investigations. 


Adaptation Liquid-Liquid Extraction 


Liquid-liquid extraction uranium from acid solutions has 
been developed and applied many countries the last ten 
years, point where would difficult review detail. 
the United States, the research groups that have pioneered this 
research, under USAEC contracts, include the Dow Chemical 
Co. California, the Union Carbide Nuclear Company Oak 
Ridge, Tennessee, and the Pont Nemours and 
Company the Savannah River Laboratory, S.C. 


The phosphate (TBP) refining process was one 
the first achieve success this field and has now attained 
high degree efficiency. Data are available from operating 
refineries predict that the process could possibly applied 
the direct treatment nitrate eluates from ion exchange 
circuits shown Figure The mechanism extraction 
favored high nitric acid concentration normal) and de- 
pressed the presence the sulphate ion. Since eluates 
produced leach plants contain essentially 0.3 0.4 HNO; 


Figure 4—Four stage mixer-settlers. 


DIALKYL ALKYLPHOSPHONATE TERTIARY AMINE 


| Recycle Borren 
Nitrate Eluote 


Barren 
Eluote 
To Leach 


Figure Proposed flowsheets for refining uranium 
from ion exchange eluates liquid-liquid extraction, 


and 0.3 0.6 these concentrations would not 
favourable equilibrium those obtained from the dissolution 
uranium concentrates nitric acid. These factors and their 
effect possible uranium loading countercurrent circuit 
would require further study. The ultimate effect the purity 
the product could best determined continuous testing 
liquid-liquid extraction apparatus. 

second possibility can also considered where TBP 
replaced dialkyl with more favorable 
extraction coefficient. Some kerosene-soluble phosphonates are 
commercially available this time and the more suitable com- 
pounds could chosen the basis low water-solubility, 
high capacity and lower cost. Available data their behavior 
were limited and little work had been reported predict 
whether water other solutions could used for stripping 
practical countercurrent system. 

third possibility appeared the adoption sulphuric acid 
uranium from anion exchange resins place 
nitrate chloride solutions. This would allow the use liquid 
anion exchangers, such the tertiary amines, for re-extraction. 
These highly selective extractants the sulphate anion 
complex are now used many leach plants place 

Ammonium diuranate (ADU) has been shown one 
the more suitable forms uranium for the fabrication high- 
density dioxide pellets. The present investigation was therefore 
concerned with the suitability the ADU which would 
produced from the proposed flowsheets Figure con- 
tinuous bench-scale circuit. 


EQUIPMENT 


Since 1948 the Mines Branch has been active the develop- 
ment leaching, ion exchange and liquid-liquid extraction 
applied Canadian uranium and thorium ores. The research 
facilities include small pilot plant section where liquid-liquid 
extraction and stripping unit was used this study. This unit 
will described briefly the following paragraphs. 

Two four-stage countercurrent mixer-settlers, shown 
Figure were constructed from Plexiglas. The extractor 
section rectangular tank 16-in. 12-in. high. Each stage 
composed mixing chamber (3-in. 4-in.) which the 
aqueous phase introduced from the previous stage through 
bottom weir and draught pipe extending from false- 
bottom plate set one inch from the bottom. The organic phase 
enters over weir from the opposite side. Mixing done with 
three-blade, 2-in. diameter marine impeller with upward thrust, 
driven variable speed motor. The impeller acts pump 
for the aqueous phase which must move against level 
differential from stage stage. The organic phase flows 
gravity from weir weir, each weir being lower than 
the preceding one. This arrangement dispenses with the use 
interstage pumps airlifts. The mixer compartment connects 
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Figure mixer-settler 


with the settler through vertical slot (4-in. and 
around L-shaped This allows internal recycling 
the organic phase that the mixing ratio may kept 
much higher value (O/A 2/1) than the actual flow ratio. 
The O/A flow ratio set the requirements the system, 
which must achieve high recoveries near saturation the 
organic phase. The extract from this section overflows directly 
stripper the same design but reduced half the mixing 
volume and half the settling area; the overall dimensions this 
case are 13-in. 6-in. high. 

Accessories these two units include solvent metering 
pump (rated 700 ml./min.) recycle the solvent, and 
control the eluate and strip solution flowrates. 
For control the strip solution, automatic recorder- 
controller and titrimeter valve were used the first stage 
mixer compartment, while the second, third and fourth stages 
were controlled manually. 

For precipitation ammonium diuranate from the strip 
solution, continuous circuit was assembled from two-litre 
mixing tank equipped with heating mantle, mechanical stirrer 
and electrodes. Metered streams solution and ammonia 
gas were fed near the bottom this tank, while the slurry 
overflowed surge tank and hence vacuum stationary 
filter. 


Procedure 
Source Materials 


Drum quantities ion exchange nitrate eluates were obtained 
request from two leach plants Elliot Lake, Ontario. The 
solvents were prepared from commercially available 
phosphate (TBP), dibutyl butylphosphonate (DBBP) and di-2 
phosphonate (DEEP), all which were 
supplied the Electric Reduction Co. Canada Limited, 
Toronto, Ontario. 

For the study the sulphuric acid elution-amine extraction 
process, the eluate was produced the Mines Branch during 
pilot plant leaching operation. Sufficient leach filtrate was 
carry out cycles adsorption and elution 
two four-inch diameter ion exchange resin columns. The sul- 
phuric acid eluate collected from this investigation was used 
feed the liquid-liquid extraction circuit. The solvent chosen 
work was tri-capryl amine (TCA) supplied General 
Mills, Inc., Kankakee, under the trade name Alamine 
336. This amine produced commercially and has had wide 
acceptance acid leach plants for example, the Eldorado 
Port Radium leach plant from 1958 1960. Previous studies 
the authors sodium hydroxide-sodium sulphate stripping 
amine extracts® had indicated certain advantages over the 
nitrate carbonate stripping procedures. For this 
was felt that, for the production ADU, ammonium 
hydroxide ammonium sulphate strip should investigated the 
present circuit. 
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Ca: ELUATE 


Figure 6—Uranium extraction equili- 
bria 25% TBP nitrate eluates. 


Ca g U,Q,/1 AQUEOUS 

Figure 7—25% TBP-uranium extrac- 
tion and stripping mixer-settler 


circuit. 


The kerosene diluent used all three flowsheets was ob- 
tained from the Shell Oil Co. Canada Limited, and the same 
refined, high-flash-point diluent used the Eldorado 
refinery Port Hope, Ontario. 


Operating Conditions 

The circuit, shown Figure was operated for periods 
six seven solvent replacement (20 litres), samples were 
taken the equilibrated aqueous phase each stage extrac- 
tion and stripping, together with samples the solvent extract 
and recycle. composite barren eluate from each run was also 
used determine soluble and entrained solvent losses. For the 
amine studies, amine determination the solvent was also 
done standard acid-base titration. 

Dilution the TBP and phosphonates 25% v/v 
kerosene was chosen the basis calculations which had 
shown that order obtain high recoveries (99*%) 
practical number stages, say four five, such concentration 
would necessary. The amine was diluted v/v and this 
dilution was dictated the problem phase separation and 
stripping higher amine concentration. The diluted amine was 
modified with n-decyl alcohol (3% v/v) further improve 
phase separation and prevent third phase formation. 


RESULTS AND OBSERVATIONS 


Equilibria and Stage Calculations 
Tri-n-butyl Phosphate 

series seven continuous runs was completed various 
nitrate and sulphate levels one feed eluate. Equilibria 
obtained from the are presented graphically 
Figure which emphasizes the counteracting effects the two 
variables the distribution coefficient (Co/Ca). 3.3 
HNOs the coefficient was depressed from increase 
the sulphate from 0.42 1.2 controlling the sulphate 
0.5 the lime precipitation step, coefficient was 
obtained with feed eluate 1.2 HNOs. 

The data curve from Figure were transferred 
Figure for graphical calculation the number stages 
required lower the uranium content acceptable value for 
recycling ion exchange barren eluates. value 0.1 
the maximum allowable for efficient elution. will seen 
that minimum five stages are needed with feed eluate 
analyzing The O/A flow ratio 1/4 was chosen 
limit the number stages practical value. order 
approach saturation the solvent (80 much 
larger number stages would needed. 

Stripping much less problem. high uranium con- 
centration the strip would preferred were the solution 
evaporated. this study, since uranium was recovered 
ammonia precipitation and the concentration was not first 
importance, the O/S flow ratio was set higher 2/3 reduce 
the number stages needed for near complete stripping. 
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Figure 8—Batch extraction equilibria 
25% v/v dialkyl alkylphosphonates 
nitrate eluates. 


Dialkyl Alkylphosphonates 


Because the limited data uranium extraction low 
nitric acid concentration, the work was directed first some 
equilibration tests separatory funnels. 25% v/v solutions 
DBBP and DEEP conditioned with nitric acid 
solution the same normality the eluate, were then equilib- 
rated with leach plant eluates two nitrate and two sulphate 
levels. The results shown Figure indicate that for equal 
concentrations volume, DBBP more effective. Saturation 
capacity double, while the extraction coefficient 
(Co/Ca) low concentration times higher. must 
noted, course, that the molar concentration solvent 
and 0.6 for DBBP and DEEP Later 
equilibration tests with synthetic eluates and 25% DBBP 
emphasized the effect the total nitrate and sulphate ion con- 
centrations constant acidity 0.4 The data are 
presented Figure which can used guide the choice 
nitrate and sulphate levels used practice. 

The batch stripping equilibria Figure were obtained 
with water alone and with ammonium sulphate complexing 
agent. Whereas water relatively poor stripping agent, 
10% w/v solution ammonium sulphate adjusted 3.5 
100 for DBBP and for DEEP. 


the basis the above results and its higher rate phase 
separation, the 25% DBBP solvent was chosen for continuous 
testing the countercurrent mixer-settler circuit. 


series nine runs three sources eluate was completed 


Figure Effect nitrate and sul- 
phate uranium distribution coeffi- 


cient using 25% v/v DBBP. 


Figure 10—Batch stripping equilibria 
25% dialkyl alkylphosphonates. 


uranium concentration obtainable was 0.9 
sulphuric acid. This solution was fed continuously the 
tion circuit series fifteen runs. Extraction and stripping 
equilibria typical run are shown Figure 12. 

this system, contrast with the preceding two, uranium 
concentration both phases much lower dictated elution 
results and extraction equilibria. From the data Figure 
can seen that the barren eluate recycled 
leaching, three extraction stages will sufficient. the 
stripping section, where pinch noticeable the stage diagram, 
the uranium recycled the extraction circuit will minimum 
the relativ ely high distribution coefficient the third fourth 
stage extraction. Also emphasized Figure the need 
for close control each stage stripping. slightly 
below the point uranium precipitation was maintained 
order take advantage the improved stripping coefficient 


the higher pH. 


Precipitation 


The strip solutions collected from the liquid-liquid extraction 
circuit were used produce ammonium diuranate (ADU). 
Batch and continuous precipitation tests with gaseous ammonia 
were done study the effects pH, temperature, and ammonia 
gas dilution filtering, washing, and subsequent reduction and 
sintering 


with 25% DBBP. Results typical run are presented graphi- 

with the TBP uranyl nitrate system, recovery was 

with w/v ammonium sulphate with the the first 

was obtained four stages. The second, third and fourth 

stages did not require control; when the stripping feed 2(pH 3-9) 

was adjusted 6.5, the these stages was slightly 

cycling 50% the barren strip, which had been done 

order the ammonium sulphate consumption. ‘The 

nitrate ion depresses the stripping and results showed 

that equivalent stripping was obtained with ammonium 

ool 10 100 


Amine 
The feed cluate was produced from the elution the pilot 
plant ion exchange columns, with 10% The maximum 
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Figure 11—25% DBBP-uranium extraction 
equilibria mixer-settler circuit. 
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Figure tri-capryl amine-uranium extraction and 
stripping operating equilibria mixer-settler circuit. 


With the TBP strip solution, where the sulphate ion 
absent, optimum conditions for good filterability and high final 
density the pellets were defined follows: 


Temperature 40-60°C. 
Retention time 20-30 min. 
Gas dilution None 


Drying temperature 


Under these conditions sound pellets with densities 10.5 
10.7 were produced. 


Batch precipitation phosphonate and amine strip solutions 
with high sulphate content produced coarse ADU containing 
15% This had been experienced leach plants 
where ammonia used precipitate uranium from the eluates. 
With continuous precipitation 7.0 7.5, this was 
avoided. the other hand, the product was very fine and 
filter and wash. This difficulty was minimized 


TABLE 


higher temperature and the dilution the ammonia with 
nitrogen. The optimum conditions were defined follows: 


Temperature 60°C. above 
Retention time 10-20 min. 


Uranium dioxide pellets, having reproducible densities 10.6 
10.7, were obtained from this ADU. 


Analyses ADU Products 


Chemical and spectrographic analyses the ADU products 
are shown Table For comparison, the analyses typical 
leach plant ion exchange product and the specifications the 
ADU product from the TBP refinery are included. terms 

purity, the leach plant product relatively 
high iron, thorium, calcium, boron, molybdenum, and rare 
earths. these, iron and calcium are easily removed any 
one the three processes. Thorium effectively reduced 
specification level, except the case TBP, which low 
degree saturation the extract will allow extraction some 
thorium with the uranium. the amine process which 
anion exchange reaction, thorium essentially rejected, but 
molybdenum, even small concentrations, extracted. Boron, 
high neutron absorber, reduced below specification. Lantha- 
nide rare earths were also reduced below specification with 
results are lacking for TBP and TCA. 

this point must stressed that the significance these 
analyses terms purity can only deter- 
mined carrying out neutron absorption tests and in-reactor 
evaluations the sintered dioxide pellets. 


Solvent Losses 

Entrained and soluble solvent losses the barren eluate 
were low per pound recovered uranium and will represent 
negligible item the overall reagent cost. With DBBP the 


CHEMICAL AND SPECTROGRAPHIC ANALYSES ADU 


Source 


Leach Plant 


Chemical 
U;0s, 70 80.0 88.0 
3.0 1.8 
NO, 0.2 
0.2 
0.5 
0.5 
Th, p.p.m 500 100 
100 
Spectrogra phic 
Ag, p.p.m. 
5-10 0.1 
<10 <10 
1750 
<10 <10 
20-100 <10 
<20 
<300 
150 
<50 <10 
220 
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TBP from 


Refinery 
DBBP from TCA from Specifications 
85.5 
2.8 
1.8 
<0.1 0.2 
1.0 
<0.1 <0.1 0.2 
0.2 
<10 <10 
<17 
<10 
200 
<10 
0.06 0.1 
0.06 0.05 
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TABLE 


COMPARISON REAGENT FOR REFINING 
URANIUM FROM ION EXCHANGE ELUATES 


Reagent 
Process 
Price Consumption 
TBP HNO; 0.07 2.5 0.175 
NH; 0.045 0.7 0.031 
Total 0.206 
DBBP HNO; 0.07 0.04 0.003 
NH; 0.045 0.01 
Total 0.0274 
TCA 0.012 2.0 
Lime 0.01 —0.015 
Total (savings —0.110 


loss was 0.006 while with TCA was too low 
determine within the 100 solvent cycles completed the 
circuit. Loss TBP was not determined. 


Reagent Consumption and Costs 


list the reagent consumptions and costs producing 
refined ADU the three processes presented Table 
These are over and above the cost producing the current 
leach plant ADU the two-stage precipitation 
nitrate eluates. the case TBP, the eluate after lime 
neutralization reacidified 1.0 and the barren 
eluate neutralized again with ammonia 0.4 before 
recycling elution. For DBBP, the only reagents involved will 
the ammonium sulphate for stripping, plus the small amount 
nitric acid extracted which must neutralized with ammonia 
during stripping. the TCA process the sulphuric acid 
recycled leaching except for the amount adsorbed the 
resin and the extract; this acid replaces the total nitric acid 
and lime consumed the present nitrate elution process and 
shown negative cost, savings. 


DISCUSSION 

the basis purity product, the dibutyl butylphos- 
phonate (DBBP) process will preferred. also more 
readily adaptable acid leach plant flowsheets, which either 
all the eluate only portion can treated without disrupt- 
ing current procedure. acidified sodium chloride elution 
use, change can readily made nitrate elution without any 
change the equipment. 

The amine (TCA) process does have some good 
points. First, the need for recycling the “royal 
uranium bleed the first fractions the barren effluent will 
eliminated when sulphuric acid elution used. Secondly, the 
elution done open circuit, which eliminates possible 


build-up impurities. Thirdly, where nitrate elution presents 
problem waste disposal because the very low permissible 
nitrate level natural streams, the sulphuric acid elution would 
easy solution. Lastly, there saving reagent costs 
but the expense longer elution cycle. noted previously, 
molybdenum the final product could major problem with 
certain ores. such case, ammonium chioride strip 
place the proposed sulphate strip would more selective for 
uranium, suggested the this would have 
confirmed for individual cases. 

The tri-n-butyl phosphate (TBP) process presents major 
increase cost and would not competitive with current 
refinery costs. The purity the product would improved 
phosphate complexing the thorium. 


CONCLUSION 


Three proposed flowsheets for refining uranium from ion 
exchange eluates have been studied sufficiently continuous 
bench scale equipment assess their feasibility. shown 
Table the dibutyl butylphosphonate extraction process will 
produce nuclear grade ammonium diuranate additional 


reagent cost three cents per pound uranium. The amine 


process affords saving eleven cents reagents but may 
not readily adaptable; will require further testing 
assure its success individual cases. The TBP process applied 
nitrate eluates not competitive with the 
applied the chemical precipitate. The work has therefore 
established that nuclear grade uranium product can produced 


the uranium leach plant relatively low cost and with little 


change the existing flowsheet. 
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Photography 
Bubble and Drop 


and 


Methods photographing bubbles and drops 
moving liquid field are described. Shape, size, 
terminal velocity, local semi-vectorial velocities, size 
distribution and high speed phenomena such 
coalescence and breakup can advantageously 
investigated photographic methods. Recommended 

lights, reflectors, diffusers, screens, 
and f-stop values are given. Applications 
particle trace and interrupted particle trace photog- 
raphy are described. 


still and motion picture photography have been extremely 
useful capturing bubble and drop phenomena form 
record which often reduces the personal factor inherent 
Such photographic records must interpreted 
with considerable reservation. For example, large air bubble 
cylindrical tube may appear exhibit ledge along 
periphery its spherical cap still photograph 
with high speed electronic flash will show such 
but disturbed, wavy periphery. Only when high 
speed motion picture camera used and the results viewed 
motion screen can the situation resolved. The 
the bubble fluttering, wavy ribbon which appears 
ledge due hysteresis the eye. The steady state ledge 
formed, but only large bubble viscous (over 
liquid. Photography has aided the resolution many 
phenomena, usually the acquisition qualitative 
The principle uses photography bubble and drop investi- 
have been record the shapes dispersed 
scale the scene, drop and bubble sizes may 
quantitatively adding dye one phase, 
the boundary layer separation can taking 
various light transmission disturbances the flow 
patterns the field fluid can studied. Such disturbances 
the photoviscosity effect bi-refringence 
fine solids such differences con- 
and changes refractive index due tem- 
Stroboscopic and other multiple image 
state and measure velocity travel. High speed 
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motion picture techniques permit observation the nature and 
progress such fast phenomena bubble drop forma- 
breakup. Much has been learned 
boiling phenomena the use such They 
also permit close views with sufficient time for observation 
without other elaborate summary readily 
available pertinent articles given Table 


Shape and Eccentricity 


Rising gas bubbles are easily 
two-dimensional shape alone desired, any camera with 
double triple extension bellows and almost any film may 
used. Backlighting best and electronic flash unit 
produce short duration light sufficient intensity give 
excellent results. good home built unit high power can 
built for fraction the cost commercially sold one. 
good images the surface the bubble are desired, 
unit top light should added. When photo flood lamps are 
used, the exposure time must the order 0.001 second. 
This means, course, either fast film very strong light 
must employed. One may assume that large gas bubble, 
film. the image the film may allowed move 
0.1 mm. and the bubble photographed full natural 
size, the time exposure but 1/3000 second. Electronic 
flash units are more economical over extended service than 
photoflood use reflectors such white 


TABLE 
PHOTOGRAPHY BUBBLE AND Drop RESEARCH 


photographs show shape, 3-9, 11-13, 16-20, 26, 


eccentricity, boundary layers, etc. 27, 29, 34, 

Size and shape shadowgraph 13, 14, 

Velocity multiple image still 

photographs 

Velocity dark field trace 29, 31-33, 

Velocity motion pictures 3-5, 14, 37, 

Drop internal velocity pattern 18, 29, 34, 

dark field trace 

Strobolume methods 10, 

Velocities interrupted trace 32, 

Shapes and qualitative action mo- 26, 27, 

General qualitative and ative 15, 17, 26, 27, 
action high speed motion pictures 36, 37, 
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Lamp 
+) 


Dif fuser 


RFL2 


Tank 


Bubble 


Light 
Shield 


Camera 


Figure Photography with back 
lighting. 


cardboard, aluminum foil, and glass mirrors can greatly aid 
good pictures with single light source. Cardboard shields 
must used protect lens from unwanted direct rays 
from the lamps from reflecting surfaces. good arrangement 
shown Figure When the dispersed particles must move 
vertical cylindrical vessel, lighting must supplied from 
such direction that images the light sources will not enter 
the camera lens reflection from the cylindrical tubes. 

Liquid drops moving liquid medium are more difficult 
photograph due the much smaller differences densities 
and refractive The simplest solution dye the 
drop produce high Calco Oil Red 1700 
excellent and seems entirely dev interfacial activity char- 
acteristics, the presence which cannot determined 
surface tension measurements alone. combination side 
and back lighting has been found best. tan light brown 
reflecting background gives good results with panchromatic 
films such Eastman Tri-X Ansco Super pan. good 
arrangement for both still and motion pictures shown 
Figure background paper must far enough away 
from the container avoid shadows. Some photographers 
prefer close background any kind. dye used the 
drop liquid, color photography can made superior black 
and white. Garner and used back lighting with 
tungsten lamp and fast shutter speed produce excellent photo- 
graphs the boundary separation and wake slow moving 
drops alcohol water. 


Drop Size Distribution 


One may photograph the dispersed droplets fluid field, 
enlarge the photograph and laboriously arrive drop size 
distribution. Care must taken avoid photographing drops 
adhering wall. The collimated flat beam illumination shown 
Figure the best insurance against erroneous results. 
this method only freely suspended drops are illuminated. 


Velocity Measurements 


Three kinds velocities must considered bubble and 
drop research. True vectorial velocities cannot photographed. 
Semi-vectorial velocities can used plot velocity fields and 
true vectorial velocities may then obtained from these 
process graphical differentiation. The most commonly 
measured velocity the “gross terminal obtained 
measuring the time necessary for particle fall (or rise) 
through measured vertical distance. The distance must 
long enough that all cyclic phenomena such 
will repeated many times. Several photographic methods 
have been developed for such determinations 


Velocity Multiple Image 


The useful device for short distance 
and time determinations. About six images can obtained 


Paper Background 


Shield 


lighting. 


Arc Lamp Black Background 


RFL2 


Camera 


Figure 3—Side lighting with 
mated beam. 


single piece film. Local variations velocity can detected 
and stages short cycle phenomena can observed. One can 
also technique whereby the camera shutter can 
operated exact time intervals cam and microswitch 
arrangement. The wear the camera excessive and this 
method not recommended. 


Particle Trace Methods 


Darkfield particle trace photography can yield much informa- 
tion local Particles carefully de-greased 
aluminum powder, clay other materials are suspended 
the continuous phase and strong flat beam light sent 
through the field right angle the line along which the 
camera aimed. timing device such motor-driven arm 
must included the photograph order measure the 
short clapsed time corresponding the streak reflected light. 
The time exposure camera shutter not 
dependable for such work. typical arrangement shown 
Figure 


Velocities Inside Drops 


The flow pattern the interior moving drop has been 
qualitativ ely recorded many times both 
have recently published very set quantitative 
measurements the velocity distribution inside falling drop, 
using motion picture technique. all these cases, the camera 
and lights were placed platform and allowed fall the 
same speed the falling drop, the rate fall the platform 
being controlled hydraulic cylinder. 

using unmarked rotameter tube with the larger end 
up, species “water was recently used carry 
limited study the circulation inside 84% 
glycerine solution was pumped upward through the rotameter 
tube. The fluid velocity was slower the top than the bottom 
and, some point the tube, its speed could made exactly 
match the fall velocity drop measured size which was 
introduced the top the system. With this apparatus the 
drop could held one point for hour more for observa- 
tion. flow diagram given Figure 


The photographic arrangement shown Figure 
from carbon arc led thru condensing lens and adjustable 
vertical slit, thus producing flat beam light illuminate 
narrow cross section the drop. square plastic tank was 
placed around the rotameter tube minimize distortion and 
keep the light beam sharp. The space between tank and tube 
was filled with 84% glycerine. The camera was set support 
90° the light beam and black cardboard background was 
provided. Degreased aluminum particles appeared the 
graphs tiny light blips and their progress recorded dark 
field particle trace photography. 
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Quick Closing 
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Assembly 
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Figure 4—Flow diagram tapered tube system. 


The camera used was Speed Graphic with triple 
extension bellows and f4.7 lens Compur shutter. Film 
used was Ansco’s Super Hypan with rated ASA speed number 
400. overdeveloping methods, its effective speed number 
was doubled 800. typical photograph shown Figure 
Several such photographs were obtained the course each 
experimental run, the total number such negatives being 
about 500. 


The problem then becomes one interpreting the streaks 
shown the photographs terms circulation velocities. 
Garner and Haycock used organic mixtures having the same 
index refraction the surrounding glycerine mixture but 
this limits the range variables involved. the present system 
light reflected from the various 
traveling along the streamlines inside the drop result 
circulation pass through the interface the drop 
fuid, through the glass the tapered tube, and finally through 
the plastic the box surrounding the tapered tube assembly. 
There difference index refraction between the drop 
fuid and the field fluid, and this difference presents the most 
error the interpretation photographs. method 
presented below. 

The assumptions involved this development are the 
following: 

(1) The curve generated the intersection plane 
passed through the vertical axis the drop can expressed 
equation form. 

(2) The drop symmetrical about its vertical axis. 

(3) The pencil rays reflected from the aluminum particle 
which impinges upon the lens the camera can approximated 
single ray ray) passing through the drop inter- 
face and moving thence direction normal the plane 
the sheet light passing through the drop. 


Figure 6—Drop not centered tube. 


Figure 5—Trace photograph drop. 


(4) The lens the camera lies along line normal the 
plane the sheet light, and which intersects the major and 
minor axes the drop. 

(5) The distortion the hardware about the drop 
negligible. 

The first assumption seems justified although its actual 
realization practice may quite difficult. 

the drop centered the tapered tube, and the tapered 
tube perfectly vertical, the second assumption justified with 
respect the photographs the drop obtained. Whether 
not the drop vertical axis lay along the tapered tube axis could 
checked this equipment suspending drop about 0.8 
cubic centimeters volume, the density which differed from 
that the field fluid about 0.1 gram per cubic centimeter, 
the tube and observing the drop’s symmetry about 
the vertical axis from two directions. drop photographed 
the tapered tube which had been thrown out adjustment for 
the purpose shown Figure This figure also shows the 
possibility photography non-symmetrical drops this 
apparatus purposefully causing the tapered tube out 
adjustment. 


The third assumption greatly simplifies the treatment the 
drop data and necessary for this reason. The error introduced 
this assumption seen quite small, the direction 
the ray considered below differs from that the 
direction the true “average” ray very small quantity. 

The fourth assumption can realized proper design 
the equipment. Note that because this assumption, the position 
the drop the tapered tube cannot allowed vary 
vertical direction. 

The fifth assumption justified the basis photograph 
spherical ball suspended the tapered tube under nonflow 
conditions. Photographs spherical balls different sizes 
taken the course the work show equally small amount 
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Figure trace photography. 
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distortion caused the optical system the tapered tube 
and the plastic box surrounding it. 

The problem one finding the actual location the 
aluminum powder trace shown the drop photograph with 
respect the geometry the system. graphical method 
available, but somewhat tedious and does not lend itself 
the processing large amounts data. more direct analytical 
method based the above assumptions also available, and 
modification this method which does lend itself the process- 
ing large amounts data discussed below. 

Snell’s law refraction familiar from the study optics. 
is: 


This must applied the plane the normal the lens 
surface (at the point intersection the ray with the surface) 
and the ray. seen that method solution could 
find the equation the surface the drop, find the normal 
the point intersection, and proceed the determination 
the deviation from the apparent position the aluminum 
particle the application trigonometry. illustration 
this shown below assuming that the curve generated the 
intersection the plane the sheet light with the drop 
interface that ellipse. 

found practice that the curve the drop interface (as 
obtained above) seems follow one equation above and another 
below its major axis. these equations are both those 


ellipses: 
(2) 
(3) 


respectively, seen that: 
a, = d& by by (4) 


for the drops this work. 
Consider the upper ellipse given Equation (2), which 
the equation given the simultaneous solution of: 


and: 


the case drops symmetrical about the vertical axis 
The apparent position the trace available from the photo- 
graph question. The actual physical dimension the major 
and minor axes the drop also available from the photograph 
and knowledge the width the tapered tube the position 
the drop. Let the length the major axis taken two 


The values and are available from the photograph. 
They must substituted into Equation(8) obtain the value 
that 


The equation plane tangent the surface the drop 
the point its intersection with the light ray from the 
particle is: 


and the direction the normal the surface the drop 
this point given by: 


2 
1 


where: 


Thus the direction the normal is: 


The angle between the line corresponding the path the 
ray from the drop interface the camera lens given by: 


arc sin 


This substituted into Equation (1) that: 


ne 
Now: 
and: 


The direction the photograph given by: 


this way several points the trace can selected and 
the actual position the points plotted the 
The length the trace can found means calipers, and 
the velocity the particle which caused the trace determined. 
Final calculation results awaits only the availability time 
digital computer complete the thousands such solutions 
necessary. 


Interrupted Trace Methods 


slotted metal disc, driven known exact speed 
synchronous motor, inserted between the 
and the subject, the slots the will enable stroboscopic 
photography the light narrow slot 
followed wide one and space provided between 
such pair and the succeeding one, dot-dash image will 
The number such images must the same for 
every sequence the particles are travelling the 
lesser number dot-dash pairs indicates that the particle 
was travelling across the beam and the data must rejected 
invalid. The known rotational speed the spinning 
disc provides the time measurement. 


Flow Around Cylinder 


large gas bubble will rise through liquid shape which 
has been described spherical The front surface 
portion perfect and the rear area disturbed 
but relatively flat somewhat concave one. very large 
liquid drop falling through viscous liquid would possess 
definitely crescent-shaped vertical cross The profile 
the frontal surface the rising bubble was determined 
Davies and and idealized geometric shape described 

The vertical cross section the idealized spherical cap will, 
when described but two dimensions, circular segmental 
cylinder. Such cylindrical shapes were made plastic and 
pulled through water rectangular plastic tank 
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Figure 8—Elevation interrupted trace apparatus. 


shown diagrammatically Figures and The light from 
arc was led through spinning disc interrupter, 
adjustable slit and condensing lens the tank. This produced 
aflat plane light permit two-dimensional field illumina- 
tion. Thus, system dark field, interrupted streak photo- 
sraphy was made possible. Much the equipment had been 

The light interrupting apparatus consisted rotating 
disc shown Figure metal driven synchronous 
motor operating 1700 r.p.m. Two slots were provided. One 
was cut out and the second, located the circumference 
away, was 30° width. The spinning disc was located front 
the adjustable slit. The exposure times are regulated the 
width the slits and the width the light beam the adjust- 
able slit. 

The camera used was Speed Graphic with 127 mm., 
f4.7 Ektar lens. The full lens opening was needed, using 
Eastman Tri-X cut film. The shutter speed the camera was 
permit least two revolutions (dot-dash images) the 
disc. The still pictures were supplemented with high 
seed (300 pps) motion pictures using Fastax WF-3 camera. 
The arc lamp was then operated direct current avoid the 
fuctuations light inherent supply. The interrupter 
was not used for motion pictures. 

The still camera shutter was automatically triggered 
means photoelectric cell and circuit. light beam was 
from mirror the bottom the weight pan 
that the shutter would opened while the rising cylindrical 
shape was the field view, shown Figure 10. 

obtain velocity values, the negatives were placed 
and projected piece graph paper. Each 
discernible sequence particle images was traced 
the paper along with the base the lowest image the shape. 
After the images had been traced, both the horizontal and vertical 


Figure trace photograph. 
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Figure 9—The interrupter disc. 


distances between successive dashes were measured 
tioned relative the center the base line the shape. 
Horizontal and vertical velocity components could then 
determined means the time elapsed during one interrup- 
tion cycle. keeping the geometry both photography and 
enlargement constant, the results considerable number 
repetitive runs could plotted single sheet graph paper. 
moving the reference center the graph paper from the 
center one shape base line that the next, the velocities 
may converted those relative the shape. originally 
photographed they are relative the tank and other stationary 
parts (Figures and 12). 

The number dot-dash images must correspond the 
number images the moving shape each photograph. 
lesser number traces was recorded, the particle was moving 
across the flat light beam and that trace was ignored. more 
than very few such invalid traces recorded, the whole 
experiment must rejected and the cause the action found. 

The method described permits the determination local 
velocities without the great expense “water tunnel”. The 
results can related the usual dimensionless groups can 
interpreted terms the hydrodynamics moving sub- 
merged bodies. 


Velocity Motion Pictures 


Davies and measured the gross terminal velocity 
rising air bubbles means periscope arrangement some- 
what like that shown Figure 13. The framing rate the 
motion picture camera measures the time interval between the 
image the bubble passes the lower mark and, later, the 
upper one. Calibration is, course, necessary. 

and Haberman and used the progres- 
sion the bubble image across succeeding frames mm. 
motion picture film determine velocity rise. 


Figure 11—Particle traces around cylindrical cap shape. 
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Figure 12—Observer moves with cylinder. 


The gross terminal velocity single drop moving 
liquid field easily any several techniques 
with error less than 4%. The velocity single drop 
moving stream such drops more usual 
method divide volumetric holdup volumetric flow 
rate obtain residence time. This then used with vertical 
distance obtain average velocity rise fall. technique 
for measuring the velocity each individual drop rain 
such drops the use ordinary (up p.p.s.) motion 
pictures 

The containing tube with its rain falling drops was 
the center the scene, 0.01 sec. electric timer the left 
and panel which record pertinent data the right 
shown Figure 14. scale was provided the rear the 
vertical tube. The whole was then calibrated for geometry 
and distortion. 


The camera used was Kodak K-100 with a.25 mm. lens. 
Eastman Tri-X Reversal, mm., film was used. Illumination 
was three 375-watt reflector floods the main subject and 
small spot the stopwatch. Any the popular photocell 
light meters and good camera stand completed the equipment 
needed. The background was tan cardboard. 

good editing machine which single frames 
observed necessity. One simply records the position 
the drop enters the view together with the corresponding 
reading the stopwatch. The drop followed visually 
film slowly fed through the viewer until second pair 
readings may taken. The gross terminal velocity then 
Ah/At. knowledge the range velocities and some 
information their distribution can obtained. 


Tank 
Drops 
Scale 
Data 
Panel 


Sec. 
Electric Timer 


Figure 14—Arrangement apparatus before camera. 


Line 


Figure 13—Velocity split field photography. 


There are some obvious deficiencies the method. One 
does not know the front-to-back position the droplet and 
velocity vs. diameter distribution tenuous. This has now 
been remedied the use side- -lighted, dark field methods 
illustrated Figure Even p.p.s. not fast enough for 
systems high density difference which the drops fall 
high velocities. The authors are fortunate having Fastax 
camera extend their range. Reflector spot are needed 
for such camera speeds. 


The results the studies which the method used are 
part series drop velocities tanks. 


High Speed Motion Pictures 


Motion picture cameras employing shutter system 
which the film stopped, picture taken and the then 
moved forward into position for the next picture are generally 
limited p.p.s. When projected p.p.s., the motion 
slowed factor four. cameras capable higher framing 
rates, the film moves through constant speed and framing 
accomplished other mechanisms. One these the rotating 
prism type, the most commonly used ones having speeds about 
200 8000 p.p.s. 

While such cameras have many uses, they are chiefly 
employed produce slow-motion movies. Such fast phenomena 
partial coalescence and can more clearly under- 
stood observation the details the action. Even 4000 p.p.s. 
was found somewhat slow observe the early stages 
bubble growth during boiling. Circulation patterns inside 
falling and the collapse steam bubbles have 
been determined such methods. The frequency and amplitude 
drop oscillations and rupture can also taken from such films. 

one wishes take motion pictures fast moving drops 
with the drop filling the field view, can use the falling 
camera, keeping the droplet view falls. telescopic 
lens and very small stop (about f-22) are necessary have 
both close view and sufficient depth field keep the 
image focus. High speed (150 500 p.p.s.) cameras permit 
one attain the same result more satisfactory 
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Diffusion Coefficients 
for the System: Ethanol-Water 


Measurements the liquid diffusion coefficients 
over the complete range compositions the 
ethanol-water system 25°C. have been carried out 
using the porous diaphragm technique with mag- 
netically-stirred and modified cell. Both the integral 
and differential coefficients agree very closely with 
the values reported Hammond and Stokes, al- 
though improved method calculaiing the differ- 
ential coefficients has been used. The limiting values 
and 100% concentrations differ somewhat from 
those previously reported. means experi- 
mental investigation the work Smith and Storrow 
the same system has been shown error. 
Analytical errors calibration have been reduced 
probable error 0.1%, although there 
average deviation for all experimental 
points relative the best line drawn through them. 
The diffusion coefficients based concentration 
driving force, show minimum mole 
ethanol, corresponding the minimum for the rela- 
tive decrease volume mixing and the maximum 
for the viscosity-composition curve. For the dimen- 
sionless Schmidt number, very sharp peaking 
value 7500 occurs the same composition. 


Ithough the diaphragm cell technique for measuring diffu- 

sion coefficients liquid systems has been widely used, 
only very few measurements binary non-electrolyte systems 
have been made. This apparent neglect has been not only due 
the difficulty constructing satisfactory cell for use with 
most organic liquids, but, more importance, the lack 
acceptable general procedure for computing the true differential 
values from the measured integral diffusion coefficients. 
obtain practically differential diffusion coefficients directly 
using small concentration differences the diffusion 
ments, extreme analytical precision necessary because the 
multiplying effects the mathematical treatment. Some doubt 
has been cast the method, also, the recent competent 
review article Johnson and Their criticism was 
based the fact that the diffusion coefficients the ethanol- 
water system obtained Smith and and 
Hammond and differ much 100%. The pairs 
investigators considered their results accurate three and 
one per cent respectively, and were not able explain the 
discrepancy between their measurements. 


Manuscript received December 8, 1960; accepted September 5, 1961. 
*Assistant Professor, Department of Chemical Engineering, Oklahoma State 
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was considered fruitful, therefore, re-investigate the 
ethanol-water system again, and possible provide modified 
diaphragm cell design and methods suitable for rigorous calcula- 
tion differential diffusion coefficients. The equations necessary 
for the latter, involving the general case for systems where 
volume changes occur have been reported else- 
The re-evaluation the diffusion for 
the ethanol-water system 25°C. presented here. 


Investigation Previous Diffusion Cells 


The diffusion cell used Smith and shown 
schematically Figure and somewhat different from the 
classical design Northrop and which provided the 
essential basis for the cells used and this work. 
Figure glass tube ending sealed male ground 
glass joint the half-cell containing the solution richer 
alcohol; the sintered glass diaphragm; are magnetic 
stirrers; large test tube serving the other half cell; 
and vented rubber stopper. 

Smith and Storrow used aqueous solutions for calibra- 
tion, determining the concentration changes the upper half-cell 
conductometrically. their experiments establish the lower 
limit satisfactory stirring rates, they found very much lower 
values r.p.m.) than (25 r.p.m.). The most likely 
explanation for this disagreement found the different 
geometry the types cells used these investigators. 
Experiments carried out the course the present work have 
established that with cells geometry similar those used 
Stokes, there are significant concentration gradients inside the 
bulk the liquids the cell-compartments stirring rates 
about r.p.m. higher rates these non-homogeneities 
cease exist. the shorter type cell-compartment, contain- 
ing the conductivity electrodes, used Smith and Storrow, 
reasonable expect that the concentration gradients the 
bulk the liquids are wiped out lower rates stirring. 
However, the electrodes, that incidentally served also 
the calibration runs, were absent the ethanol-water experi- 
ments which might have resulted inadequate stirring this 
latter case. 

The seal, essentially ground glass stopper, used Smith 
and Storrow for the upper half cell was also investi- 
using, the first instance, the simple apparatus 
Figure U-tube plugged one end with ground glass 
stopper and fitted with piece narrow capillary the 
other end prevent both eva and the development 

was filled with over 90% alcohol, was inserted, the 
edges and all the interior were dried, 
amounts water were pipetted into the bottom was 
placed The meniscus the right arm was 
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Figure 1—Sketch diaphragm diffu- 
sion cell used Smith and Storrow. 


observed move downwards while the meniscus was 
moving upwards. This process was observed for several days, 
when the liquid orignially pure water, had increased its 
volume several times the original value. This experiment 
was repeated several times although occasionally the process 
observed would not start, apparently stopped after proceeding 
far. The ground glass stopper was always inserted very 
firmly prevent ordinary leakage through the joint. The 
mechanism this distillation effect has not been elucidated 
further, although was observed that whenever the distillation 
Was progress, the edges the joint, originally dry, were seen 
moist. 

attempt measure the magnitude the distillation 
effect under conditions similar those used Smith and 
Storrow required cell, Figure that differed from the one 
used the above authors only that the sintered glass disk 
has been replaced short piece narrow capillary This 
permitted bulk flow from into case part the liquid 
were removed from via distillation through Other things 
being equal, the diffusional transport through about one 
millionth that through sintered glass ‘disk and hence 
negligible. This cell was accommodated, along with 
another reference-cell the same construction, but not contain- 
ing any ground glass joint constant temperature bath 


was filled with alcohol solution, inserted, the outside 


the cell dried carefully, about ml. water placed and 
placed position. The reference cell was treated the same 

Runs hours duration were carried out, and 
their conclusion samples were taken from using hypodermic 
needle that was inserted through and served the purpose 
venting throughout the runs. Analysis refractive index 
was carried out. The reference cell containing ground glass 
joint exhibited negligible transference, but significant changes 
were found occur the other cell and were judged due 
distillation through the joint. 
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Figure 2—Apparatus used investi- 
gate distillation through ground glass 


Figure 3—Apparatus used pseudo 
diffusion experiments. 


Several runs the low concentration differences used 
Smith and Storrow® their experiments showed that the dis- 
tillation effect might have accounted for about and 15% 
respectively, the values they obtained for the diffusivity. 
This indicates strong possibility that the higher values the 
diffusivities obtained these authors are consequence 
mass transfer through the ground glass joint used their 
apparatus. 


Apparatus and Procedures 
Equipment 

The diaphragm cell designed and used the present 
experiments shown Figure this design the use 
rubber stoppers and lubricated unusually well-fitting ground 
glass parts has been eliminated. filled and emptied through 
two narrow (0.5 mm. i.d.) capillaries, and fitted near the 
ends with valves” made from short pieces poly- 
thene tubing and commercial screwclips. These valves showed 
detectable leakage water under vacuum about mm. 
Hg. sintered glass diaphragm. are iron-in- 
polythene stirrers. the top half-cell sealed with Teflon 
stopcock. The two capillary side-arms and are situated 
the same side the cell although for the sake clarity 
Figure shows them opposite sides and set further from the 
cell- body than actuality. About the level the side-arms 
are reinforced with few windings plastic insulating tape. 

multiple unit was designed accommodate six cells 
arranged with their centres the corners regular hexagon 
(Figure 5). 

The body each cell surrounded brass sleeve 
containing four slots for aiding circulation the bath liquid. 
The cells are supported piece cork placed the sealed 
bottom the brass tube, and cork wedges the sides. 
the outside the bottom the brass sleeve there are two 
pegs, the one the centre little longer than the other one. 

The sleeves are held position bottom supports 
consisting short pieces brass rods soldered brass base 
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Figure 4—Diaphragm 


plate hexagonal arrangement. are two holes 
through each support matching the two pegs the bottom 
the sleeves hol ding the cells. 

All stirring achieved pair rotating magnets 
arranged around each cell with the poles level with the dia- 
phragms. The magnets are standing upright mounted pinions 
driven central drive and rotating bearings prov ided 
laminated phenolformaldehyde resin) sheet which supported 
centrally pillar consisting large brass tube soldered 
the base plate the this pillar there foot- 
bearing which the axle the centre gear rotates. The 
latter driven variable speed motor. The battery six 
cells assembled outside the bath shown Figure 

The constant temperature bath containing thin paraffin oil, 
was maintained 25°C. measured with standardized thermom- 
eter, and the temperature was found constant within +0.01°C. 
inside the cells. 


Diffusion Experiments 

The cells were calibrated using 0.1N. aqueous solutions 
made with Reagent Special grade and high 
quality distilled water. (Specific conductivity mho. 
the solutions were carefully degassed prior 
introduction into the cells. The solutions were degassed 
boiling, and the aqueous ethanol mixtures freezing and 
subsequent thawing under vacuum. The solutions were intro- 
duced into via (Figure was removed from the 
sintered glass diaphragms boiling the solutions under 
vacuum. Subsequently the cell was placed the bath, and after 
temperature equilibration the other solution was introduced into 
After two three hours the solution was pipetted out, 


Figure 5—Cell support and stirring device section. 


was rinsed and finally filled with fresh solucion from the 
same stock. the case the calibration runs the liquid 
initially was always distilled water. The liquid displaced 
the seal around the ground glass plug served prevent evapora- 
tion the liquid through the joint. During the diffusion 
period, stirring was carried out high rates (about 100 
After several days diffusion the contents were pipetted 
out, and the contents were blown out through 
applying dried compressed air The first few mls. were 
discarded, and the remainder kept for analysis. 


Analyses 


The aqueous solutions were analyzed modification 
the procedure recommended McBain and i.e., 
evaporating the solutions dryness and weighing the residues. 
careful operation the accuracy (+0.02 mg.) the semi- 
micro analytical balance was fully exploited. The 
change buoyancy the weighing bottles from the time 
the empty bottles the time eighing the residues 
was noted weighing “standard each time. The 
true weight the potassium chloride residue, was computed 
the following formula: 


B+R _ Me 
where 
apparent weight the bottle plus residue, 
apparent weight the bottle, 


weight standard bottle corrected vacuo, 
apparent weight standard bottle time weigh- 


ing the residue, 
apparent weight standard bottle time weigh- 
ing empty bottle. 


Three parallel determinations sample were made. 
The probable error the determinations was found 
the precision claimed McBain. 

The aqueous alcohol solutions were analyzed density 
using ml. Sprengel pycnometers modified drawing out 
the tip one the arms capillary about 0.1 mm. 
and part the other arm less than mm. i.d. the latter 
section mark was engraved. density was 
determined weighing previously calibrated standard tare 
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Figure six diaphragm cells. 


along with each pycnometer. This standard tare was prepared 
sealing off the neck 100 ml. flask. Each sample was 
subjected least two density determinations. The probable 
error the individual determinations was estimated 
g./ml. The accepted density data the U.S 
National Bureau were used. 

The commercial absolute ethanol used was treated with 
activated charcoal and fractionated reflux ratio 10:1 
ft. glass column, packed with glass helices. Its purity was 
checked density (0.79140 20°C. for 99. 35% eight) 
accurately preparing series dilutions comparing 
their compositions found the basis the measured densities 
with the calculated ones. 


The volumes the cell-compartments and the pore volumes 
the diaphragms were determined volumetrically, using de- 
gassed distilled water. 


Results and Discussion 


The analytical results for two calibration runs, Runs and 
using aqueous solutions are given Table for five 
cells, and one confirming result for Cell 


The cell-compartment volumes are given Table 2*. Each 


value the mean least three determinations. The pore 
volume for the diaphragm was found approximately 0.3 ml. 
cases. 

cell constants, were calculated the well-known 
simple logarithmic formula 


No. 6908 with the ADI Auxiliary Publications Project, Photoduplication 
Service, Library of Congress, Washington 25, D.C. A copy may be 
obtaincd by citing the Document No. and by remitting $1.25 for 
Photoprints, or $1.25 for 35 mm. microfilm. Advance payment is 
Make cheques money orders payable to: Chief, Photo- 
duplication Service, Library of Congress. 
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Figure 7—Integral diffusivity vs. mean concentration for 
ethanol-water system. 
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Ne 


100 mi. 


Figure 8—True diffusivity vs. concentration for ethanol- 
water system. 


The values recommended were used. The 
results are given Table 3*. 

The ethanol concentrations obtained the runs with the 
ethanol-water mixtures 25°C. are shown Table The 
densities used calculate the concentration values were the 
mean least two parallel determinations (exceptions: runs 
and For the column the first number 
refers the cell, the second the run chronological order. 
Two experiments with the same number represent parallel 
experiments. apparent from the table, the solution richer 
alcohol was contained the upper half-cell the majority 
runs. Owing the volume changes taking place the open 
compartment the cell the initial concentration the lower 
compartment had computed two material balances: 
one overall balance, and one either component. The necessity 
for trial-and-error solution has been eliminated simple 
device described With the exception runs 5/1, 
4/4, 1/6 and 6/7 the upper half-cell was open all diffusion 
experiments. 

From the data Table the integral diffusivities were 
computed Equation (2). Although Equation (2) has been 
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derived the assumption constant solution volumes has 
been shown that under the conditions the 
present experiments the correction necessitated the volume 
changes taking place diffusion negligible. The present 
measurements were found excellent agreement with 
those Hammond and The two sets data are 
plotted Figure integral diffusion coefficients against the 
mean concentration. 

obtain the true differential values, Hammond and Stokes 
used the method least squares after dividing the whole con- 
centration range into three regions and making rather arbitrary 
assumptions concerning the degree the power series used 
describe the individual regions. 


appeared much more straightforward use graphical 
method very similar the one used earlier which 
thought unusable the present case owing “the great 
change the diffusion coefficient with concentration, together 
with the lack any theoretical limiting values the end the 
concentration obvious, however, that the method 
least squares does not supply the theoretical limiting values 
the diffusivities but only constitutes particular way 
extrapolating the limiting concentrations. 


Plotting the combined data obtained the present work and 
those Hammond and Stokes against the mole fraction showed 
linear relationship near the two ends the concentration range. 
Hence, the limiting values could obtained linear 
extrapolation this graph. These additional points were used 
construct the best-fitting curve through the measured integral 
diffusivities shown Figure Integral diffusivities were 
calculated tabular integration for each run listed Table 
using the well-known formula 


where the bar over the concentrations denotes the mean the 
initial and final values. Next, the concentrations corresponding 
the calculated integral diffusivities were read from Figure 
and the measured integral diffusivities were replotted against 
the new concentrations Figure give the true differential 
diffusivity-concentration relationship. The very small change 
compared with Figure clearly indicates that second 
approximation was not necessary and also refutes the argument 
Hammond and Stokes concerning the difficulty using the 
graphical procedure because the very change diffu- 
sivity with concentration. any case, not the rate change 
the diffusivity with concentration, but the degree curvature 


the vs. curve that can increase the labour involved the 
approximations. the curve linear, the measured integral 
diffusivities are rigorously equal the differential values the 
mean concentrations the experiments. far the uncertainty 
introduced the use the extrapolated limiting values con- 
cerned should noted that, fortunately, their effect the 
rest the values obtained the graphical procedure decreases 
rapidly one moves away from the limiting compositions. 

Partly owing the difference the computational method 
and partly due the larger number data considered the 
present work, the resulting vs. relationship somewhat 
different from the one obtained Hammond and Stokes. The 
diffusivity values interpolated the curve shown Figure 
have been tabulated and are shown together with the correspond- 
ing figures obtained Hammond and Stokes Table 5*. 

noted that the limiting values reported this 
work 1.220 each end the concentration range compare 
the Hammond and Stokes values 1.240 and 1.132 for 
and 100% ethanol respectively. The difference due the 
additional amount data available and the improved method 
extrapolation described above. result considered that 
the end values 1.220 and 1.220 represent the best available. 


footnote for Table 
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view the close agreement between the present results 
and the data obtained Hammond and Stokes, the one hand, 
and the very significant disagreement between the above two 
sources and the measurements Smith and Storrow, the 
other, seems most likely that there was some systematic 
source error present the experiments the latter re- 
searchers. The fact that they carried out experiments 
various temperatures without any apparent intersection the 
resulting lines gives support the contention that their source 
error must have been systematic. One such possible source 
error that distillation can occur through wetted ground 
glass joint discussed above. 

Careful and rigorous examination the sources error 
and the effects the calculated results have been made. has 
been found that the most probable error single diffusivity 
value, based the analytical error, about has 
also been found, however, that the average deviation the 
experimental points from the best line drawn through them (see 
Figure Inthe case the calibration runs, the average 
deviation from the mean the individual determinations the 
cell constant equal the most probable error calculated 
from the analytical error 

There are several possible sources error that might account 
for the increased scatter the alcohol experiments. One 
them undoubtedly the fact that the thermal expansion 
cient the alcohol solutions considerably larger than that 
pure water. Consequently, equal irregularities bath tempera- 
tures give rise greater bulk flow through the porous diaphragm, 
this way the other, with the ethanol solutions. Another 
possible source error relative loss ethanol sampling 
owing its high volatility. further source error the 
fact that the diaphragms were not perfectly horizontal position 
which, under the effect larger density difference, could have 
caused some flow through the diaphragm. All these sources 
error could reduced making refinements the experimental 
technique. 

appears that these sources error played less important 
role the experiments Hammond and Stokes. The scatter 
their data seems determined primarily the analytical 
errors. This may concluded from the fact that their parallel 
determinations agree much better the range low and medium 
alcohol concentrations, where the experiments were done under 
conditions considerably greater concentration differences than 
the case the more concentrated solutions. 

Since, the present work, there was detectable trend 
the results re-runs, safe conclude that the values 
the cell constants did not decrease significantly the course 
the experiments. Also, Cell No. had been used about 
dozen runs before completing the calibration with the 
solutions without any noticeable trend the values obtained 
for its cell constant. 

Inspection the figures will show that the minimum the 
ing three molecules water and one molecule ethanol, 
very nearly so. The exact alcohol concentration corresponding 
this molal ratio 0.423 g./ml. 25°C. 

some interest note that the relative decrease 
volume mixing computed from the very accurate density 
data, found that this function also has minimum the 
same composition. 

The viscosity-composition function also has maximum 
this composition and the shape the curve reciprocal viscosity 
vs. composition also similar the corresponding diffu- 
sivity curve. notable difference, however, that the viscosity 
pure ethanol some 25% higher than that pure water. 

The fact that the above three functions have maximum 
minimum composition corresponding three molecules 
water and one molecule ethanol suggests that this composi- 
tion the ethanol-water system state maximum stability 
minimum energy. Whether not such complexes also exist 
compositions different from 0.25 question that should 
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methods investigation more suited this type 
problem than diffusivity measurements. the same time, 

when dD/dC plotted against the resulting plot consists 
two branches intersecting the composition corresponding 
the hypothetical complex, indicating that the diffusivity 
function different degree the ethanol concentration 
the two sides the minimum. Furthermore, when the diffu- 
replotted against the mole fraction the hypothetical 
complex, the plot consists two almost perfectly symmetrical 
halves the two sides the minimum, with the diffusivity 
changing only few percent from pure complex about 0.5 
mole fraction complex both directions. seems reasonably 
certain, however, that when approaching pure water pure 
ethanol other modes association become more stable, resulting 
rapid change the value diffusivity. The number 
different sorts possible associations the ethanol-water 
system being enormously large, explanation seems possible, 
the present time, the peculiar fact that the values the 
diffusivity the two ends the concentration range are practi- 
cally identical. 

Further interest terms using these diffusivity data 
mass transfer calculations arises the calculation the Schmidt 
number (dimensionless group, for the ethanol system. 
Using the diffusivity values Table with accepted data for 
viscosity and density, 25°C., found that the Schmidt 
number has values 736 and 1144 for the and 100% con- 
centrations respectively, but shows remarkable sharp peak 
about 7500 mole percent ethanol. result 
obtained the activity-based diffusion coefficient used. This 
defined (4) as: 


O In a4 
(4) 


and for mass transfer purposes highly preferable since 
based chemical potential the driving force rather than 
concentration. Calculation prov ides the same terminal 
values, and the same sharp peaking occurs mole percent, 
this time value approximately 2200. 
previously, the postulate complex tentatively advanced, 
but further elucidation terms molecular theory will 
required. 
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Nomenclature 


area diffusion, cm.? 


concentration, gms./ml. (or units shown) 

differential true diffusion coefficient diffusivity, 
cm.? 

activity-based diffusion coefficient defined Equation (4) 

integral diffusion coefficient, cm.? 

length diffusion path diaphragm, cm. 

weight, gms. 

weight standard item, gms. 

volume cell compartment, c.c.’s 

cell constant 

time, seconds 

viscosity, poises 

density, gms./c.c. 


o,f (subscript) quantities evaluated initial and final conditions 


(superscript) quantity pertaining closed and open half- 
cell, respectively 
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Concurrent Two-Phase Upward Flow 


Air and Water Through Open 
Vertical Tube and Through Annulus' 


The concurrent vertical upward flow air and 
water has been made various rates and ratios 
lus. Pressure drops and in-place ratios gas liquid 
have also been observed order correlate such 
data with corresponding flow types and patterns. 
Motion pictures taken normal and especially those 
high-speed have made possible see flow 
mechanisms indistinguishable the naked eye. 


are many industrial operations which gases and 
liquids flow vertically upward through open tubes well 
Although these operations have taken place for more 
than century, has been only within the past two decades 
that extensive laboratory studies have been made iron pipe 
well clear plastic and glass 


The visual studies which have been carried out tubing 
ranging from 2-in. size have thrown much light 
the various flow types which probably take place deep 
wells, well plant processing equipment. 

search the literature revealed the fact that little, any, 
visual studies had been made annuli, hence the decision 
compare the performance annulus with comparable flow 
rates and gas-liquid ratios open tube along with the added 
feature taking slow-motion movies both flow channels. 


Description Apparatus 


Figure can seen that the apparatus somewhat 
similar the type flow assembly which has been used 
previous investigators with provisions take such measure- 
ments as: 


(1) Flow rates both gas and liquid phases. 
(2) Shut-in ratios liquid gas. 
(3) Pressure drops over the flow section. 


Discussion Experimental Results 
Open 2-in. Tube 


Figure seen that the pressure drop decreases 
the mass velocity the air increases rate 


1Manuscript received February 15; accepted July 12, 1961. 
Petroleum Company, Borger, Texas, U.S.A. 

Chemical Engineering, College Engineering, University 
Oklahoma, Norman, Oklahoma, U.S.A. 

Based paper presented the National Meeting, American Institute 
of Chemical Engineers, Tulsa, Oklahoma, September, 1960. 

Contribution from the College of Engineering, University of Oklahoma, 
Norman, Oklahoma, U.S.A 


open tubing for low water rate For 
higher liquid rates, the pressure drops also decrease minima 
slightly higher air rates. the left these minima the flow 
type primarily that slugging. very low air rates 
not shown the graph, bubbles with occasionai bullet-shaped 
pistons work their way slowly upward through the liquid. 
this region extremely low air flow rates, was difficult 
maintain steady-state conditions from the standpoint pressure 
drop well flow rates air and water. With air 
water flow the pressure drop would, course, the liquid 
head water standing the pipe, 62.4 liquid 
height, with air only 0.076 

the right the minima, spray flow sets with frothy 
annular rings air-water bubbles and droplets moving 
relatively slow rates along the walls the pipe. the central 
spray water droplets, the velocity was much faster. These 
visual observations were clearly perceptible only through the 
showing slow-motion movies resulting from high- -photography 
800 sec. When shown the screen, the motion was 
slowed down 24/800 1/33 normal motion. 


Flow through the 5/16-in. Annulus 


Figure easy see that the pressure drop curves 
sweep through minima, however, lower air mass velocities 
than those encountered the open tube with slugging 
the left and spray-annular flow the right. With larger 
wetted-wall area per unit area cross-sectional flow, one would 
expect lesser tendency for slugging annuli. 
observed, however, that the minima for pressure drop drift 


toward the left with increasing water rates rather than the 


WATER INLET 


OPEN TUBE TEST SECTION 
ANNULAR TEST SECTION 


AIR INLET 


ORIFICE METER TRAP 


METER 


Figure 1—Flow diagram test apparatus. 
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Figure 2—Flow air and water through the open 
tube. 


right was the case with the open tube. satisfactory 
explanation can given yet for this reversal trend. 

The higher pressure drops the annulus over those the 
open tube with the same cross-sectional area can accounted 
for from the mean hydraulic radius, namely the cross-sectional 
area divided the total wetted-wall perimeter the inner and 
outer tubes. 


and Studies 


the concurrent flow gas and liquid through pipes 
various positions, has been found that the gas generally flows 
higher linear velocity than the liquid. This fact can 
proved “experimentally shutting the flow section during 
run and measuring the amount each phase place situ. 
Figures and graphs present comparison the place 
the flowing ratios liquid gas weight. 


engineers make these comparisons order 
obtain slip velocity the velocity gas compared that 
the liquid. Other investigators have termed these flow phenom- 
ena sweep efficiency curves stating that the fractional 
sweep efficiency the flowing over the place ratio liquid 
gas. For example Figure the sweep efficiency for water 
flowing has fractional sweep efficiency 
0.05 when the flow ratio 2.5 and the corresponding place 
ratio higher liquid flow rates, the sweep efficiency 
considerably greater for the same flowing ratio. 


the case the annulus (Figure 5), the sweep efficiency 
considerably higher, especially the higher water rate 
where the two ratios approach each other 
the no-slippage line equal linear velocities gas and liquid. 


Efficiencies 
The engineer primarily interested the conservation 


energy all industrial operations. For this reason Figures 
and were prepared order show the lift efficiency the 


= 0 Las/sec-rr* 


TOTAL PRESSURE DROP PER UNIT 


air and water through the 
annulus. 


RATIO 


RATIO 


Figure 4—Flow air and water through the open 
tube. 
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Figure 6—Flow air and water through the open 
tube. 


ratio the actual work lifting the liquid vertically compared 
the work done the air expanding isothermally from the 
base the top the pipe section. 


Over the range ratios and flow rates which were studied 
this experimental investigation, the efficiencies proved 
higher the open tube than they were the annulus. These 
varied inversely with the gas-liquid ratios these 
ranges; however, earlier study“ showed that the lift efficiency 
curve went through maximum 0.42 low ratio gas 
liquid and then fell off zero with flow gas. 

Strange may seem, the lift efficiency runs inversely 
proportional the sweep efficiency values. other words the 
lift efficiency was high when slugging was taking place and the 
making numerous trips and down the tube; 
however, the air was being recompressed intermittently the 
downward slippage the water giving the air more than one 
chance push upward against the water. 


the spray-annular side the picture, better sweep-outs 
were being obtained; however, much the air was flowing 
upward between the droplets the spray without doing any 
actual lifting the water, hence the lower lift efficiencies with 
the higher gas-liquid flowing ratios. 


Motion-picture Studies 

One the principal objectives these vertical flow studies 
has been that obtaining clear picture the flow types 
patterns which take place the tube annulus with various 
flow rates and ratios gas liquid. Unfortunately the written 
word even series photographs cannot convey 
adequately the visual story the reader the published page. 
Neither can the naked eye see what taking place when rapid 
changes are taking place under certain flow regimes. For this 
reason during each flow, second camera 
took concurrent motion picture record the flow 800 
frames/sec., thereby stretching out time factor when 
displayed the viewer the normal rate frames/sec. 
High- speed film and intensive lighting were required order 
give the required visibility the movie 


Description Flow Types Open Tube 


For the benefit the reader who will not have the oppor- 
tunity the film, attempt will made present 
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word picture the flow types which took place open 
1.D. tube, water was moving vertically upward 
open tubing. Air was first passed through 
the tube with the water, very slow rate, appearing 
dispersed bubbles with occasional bullet-shaped piston, the 
round nose being the top the slow-moving piston. normal 
motion view with the naked eye gives one adequate concept 
this flow type. 

Holding the water rate constant the air 
rate was increased 0.5 causing pronounced 
slugging shown Figure the top half the tube 
left-hand side the photograph, the upward moving frothy 
slug can seen, whereas the water the bottom half the 
section slipping downward along the wall the tube soon 
meet frothy slug moving upward the center the tube 
where the linear velocity higher. Thus one can observe under 
slow motion, downward slippage water along the inner tube 
wall and simultaneously upward movement the frothy mixture 
air and water the center the tube. The one-foot white 
strip has been placed the right the tube, make 
possible time the period for the slug move one foot under 
slow motion. 

With the same water rate (20 the air rate 
again increased causing spray-annular flow 
shown Figure 8B. The water droplets the spray are 
moving upward estimated velocity 100 ft./sec. through 
the center the tube, whereas small fraction the total 
water stream moving along the inner wall the pipe 
annular frothy rings linear velocities, ranging from 5-10 


Flow Types Annular Flow 


Let maintain the same mass velocity water the 
annulus, namely ft.2, did the open tube. 
air upward very slow rate, dispersed bubbles 
work their way upward through the water but the central inner 
tube (plugged both ends) prevents the formation bullet- 
shaped pistons. Figure slugging takes place 
2.5 for the air flow rate, with frothy 
slugs moving upw the center the annulus, while down- 
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Figure spray-annular 


Air rate annulus. Water rate 


strip one foot height. 


ward slippage water occurs simultaneously along the wetted 
walls the inner and outer tubes, especially the trailing end 
each slug. 

With the same fixed-mass velocity water, the air rate 
increased 16.67 with the result that spray- 
annular flow takes place did with high air rate the open 
tube (Figure 8C). 
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Analysis the Radiant Heat Absorption 
the Boundary Layer Surrounding 
Evaporating Drop’ 


approximate analysis made the effect 
the absorption thermal radiation the temper- 
ature gradients surrounding evaporating droplet 
which suspended high temperature environ- 
ment. The Hottel-Cohen method for calculating the 
radiant heat flux distribution enclosure con- 
taining absorbing medium non-uniform temper- 
ature, has been applied the spherical symmetrical 
model the system. Conduction and radial move- 
ment the cold vapor from the droplet have been 
included. 


The analysis shows that, atmospheric pressure, 
the absorption radiant energy the boundary 
layer very small and hence does not affect the 
temperature gradient nor the rate heat transfer 
the droplet surface. 


This method analysis, with the associated solu- 
tion the equations describing the reception factors, 
quite general and can applied any spherical 
enclosure containing absorbing medium whether 
one concentric hollow spheres. 


previous the evaporating rate stationary 
liquid droplet suspended high temperature atmosphere 
was shown affected the evolution the cold vapor. 
additional heat load the boundary layer was imposed 
the sensible heat requirements the evolved vapor and, 
result, the temperature profile was changed from what would 
expected there were mass transfer. This effect most 
important under free convection conditions, since the flow pattern 
around the drop dictated largely the buoyancy forces, 
which turn arise virtue the different densities tem- 
peratures the boundary layer. Any change temperature 
distribution therefore expected change the evaporation rate. 

Since many cases the vapors the boundary layer 
evaporating drop are partially opaque thermal radiation, 
might advanced that the absorption radiant energy 
these vapors might change the temperature profile the boundary 
layer and thus affect the evaporation rate. 


has shown that the absorption radiant heat 
nondiathermanous (or heteropolar) gas will seriously affect 


received January 24; accepted July 25, 1961. 
*Department of Chemical Engineering, McMaster University, Hamilton, 


nt. 
*Pulp and Paper Research Institute Canada and McGill University, 
Montreal, Que. 
Contribution from the Pulp and Paper Research Institute of Canada and 
McGill University, Montreal, Que. 


252 


Les auteurs recherchent 
tion rayonnante sur les gradients tém- 
pérature voisinage d’une gouttelette qui 
dans milieu température élevée. applique 
systéme méthode Hottel-Cohen pour calcul 
distribution flux d’énergie rayonnante dans une 
enceinte contenant milieu absorbant tempéra- 
ture non uniforme. tenu compte condition 
vapeur froide ainsi que son déplacement 
radial partir gouttelette. 


montre qu’a pression atmosphérique 
limite est trés faible par suite n’affecte pas 


Cette méthode d’analyse, qui inclut solution des 
équations décrivant les facteurs réception, est 
assez générale peut s’appliquer toute enceinte 
sphérique contenant milieu absorbant, qu’il soit 
formé sphéres concentriques vides. 


the heat transfer processes, especially high temperatures and 
pressures. Under certain conditions, heat transfer rates may 
change much more rapidly than the temperature gradients. 
recognized that absorption radiation the boundary 
layer may affect the skin temperature bodies high speed 
flow. proceeded solve the Fourier-Poisson conduction 
equation, with radiation absorption term included, but found, 
however, that the effect was small under his conditions. 

the purpose the present analysis determine the 
magnitude and distribution the radiant energy absorbed 
the boundary layer surrounding evaporating droplet. 


Analysis Problem 


formulating the equations and boundary conditions 
describing the flow and temperature distribution around 
absorption radiant energy the boundary layer cannot 
assumed negligible. The complete energy equation 
can thus written spherical co-ordinates: 


where the heat absorbed per unit gas volume. Equation 
(1) implies steady-state conditions. Because its non-linearity 


The Canadian Journal Chemical Engineering, December, 1961 


and that 
usual pre 
spherical 
and radi 
leads 
any 

Gold 
Equation 


the solut 
consider 

The 
assigned 
vary 
these 
can 
may 


(1/ra 


where 


lated 

amount 
but 
the 

bounda 
conside 
radiatic 

radiant 


Bri 
numbe 
the 
each 
can 
the rac 
Hottel 
due all 
gas. 


*The 
heat 
500°C 
is sati 


*°The 
radia 
chan 


The 


& 


and that the other boundary layer equations, direct solution 
impossible and certain simplifications must made. The 
usual procedure neglect convection effects; this implies 
spherical symmetrical model which allows radial flow vapor 
and radial temperature gradients only. The spherical model 
leads relatively simple solution, which can used justify 
any assumptions concerning the absorption radiation the 
boundary layer, posteriori. 


and have shown that this model allows 
Equation (1) written the integral form: 


the solution which more convenient, particularly when the 
gas properties are not taken mean film conditions but are 
considered change with temperature. 

The simplest solution obtained when the absorption 
radiation the gas neglected, the heat capacity the gas 
assigned average value and the thermal conductivity assumed 
these conditions, has already been that Equation (2) 
can integrated from the drop radius any point (which 
may the effective radius the boundary lay yield: 


Equation (3) allows the temperature distribution calcu- 
lated for any given set conditions and evaporation rate. 


calculate the effect absorption radiant heat the 


boundary layer, not enough know only the overall 
amount heat which transferred the gas radiation, 
but necessary know also the distribution this absorption 
inthe boundary layer. The net transfer heat any spherical 
gas shell surrounding the drop will depend not only the 
temperature this shell, but also its volume. 
Since the temperature varies continuously throughout the 
boundary layer, gas shells differential thickness should 
considered. However, because the integral character 
finite thickness must considered. 
The method solution adapted here extension the 
method presented Hottel and for calculating the 
radiant heat flux and distribution enclosure which contains 
temperature-varying absorbing gas. 
Briefly this method consists dividing the system into 
surface-gas zones, the size hich depends primarily 
the temperature variation throughout the system. Although 
zone should isothermal, some variation temperature 
can tolerated long representative temperature can 
assigned it. The procedure then write heat and material 
balances for each zone turn. This requires knowledge 
the radiant interchange view-factors reception factors (using 
Hottel and Cohen’s terminology), for each the zones, making 
due allowance for the absorption radiation the intervening 


*The latter two assumptions are reasonable since for water vapor, the 
heat capacity varies from 0.455 at 100°C. to 0510 (cal.)/(g.)(°C.) at 
500°C. and the thermal conductivity over this same temperature interval 
is satisfactorily represented by the Equation k = (57 x 10-®) + (19.5 x 


heat transferred radiation any one point the sum all 
radiation to that point from all radiators in the system. Therefore, any 
change at one point influences the heat transfer in the whole system. 
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Specifically, the problem under consideration involves the 
container which very large compared the size the droplet. 
Moreover, the temperature the wall the same that the 
gas. Under these conditions, the following two facts are 
immediately apparent: 


(1) All radiation absorbed the gas outside the boundary 
layer must re-emitted, that attenuation radiation 
results this medium. consequence this, the radiation 
reaching the drop and its immediate environment the same 
that transmitted this region were suspended non- 
absorbing 


(2) since the area the container walls very much 
larger than that the droplet, the radiation reaching the droplet 
independent the wall emissivity. 


Hence the droplet and the boundary layer only 
black environment. The model now takes the form inner 
droplet surrounded concentric sphere with radius equal 
and with absorptivity (emissivity) equal one. The outer 
surface completely pervious the radial flow vapors from 
the evaporating drop. 


The intervening gas volume between the drop and the outer 
imaginary surface, which constitutes the boundary layer the 
drop, was divided into number spherical gas shells according 
subscript number equal its radius relative the radius the 
drop. Similarly any gas zone indicated the subscript 
number equal its inner radius relative the droplet radius. 
Thus gas zone the spherical shell volume contained between 
equal and 2.5. This nomenclature shown Figure 
The inner volumes are smaller thickness, since the tempera- 
ture variations are most extreme these regions. 


formulating the reception factors for the system this 
method solution, the key parameter the absorption coefh- 
cient the gas, which measure the absorption strength 
the absorbing medium. For grey gas, constant over the 
full wave-length spectrum and the expression relating the emis- 
sivity and absorptivity gas the absorption coefficient 


where the thickness the gas. 


Absorption and emission thermal radiation real gas 
occur only over part the wave-length spectrum. 
has thus proposed that actual real gas can replaced 
grey and clear gas. The grey gas, with absorption coefficient 
can thought occupying fraction the spectrum, 
and the clear gas occupy ing the fraction a). 
essence, this scheme one fitting emissivity-versus-path 
length curve exponential equation the form: 


The parameters and depend the conditions and dimen- 
sions the system under consideration. most industrial 
systems hich the surface reflectivity low, these variables 
can found, fitting Equation (5) the gas emissivity data 
equal the mean beam length,* and such that: 


has shown that the mean beam length for any enclosure 
approaches zero. the system under considera- 


mean beam length the radius hemisphere gas which 
absorbs the same radiant flux from a surface located at the centre of 
its base, as the gas in the enclosure absorbs from the radiating surfaces. 
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Figure 1—Diagram showing spherical 
surfaces and gas zones upon which 
reception factors are based. 


tion the product very small and hence 0.667 D,, 
where the diameter the effective boundary 
Hence Equation (6) becomes 


The effective diameter the boundary layer can calculated 
from Equation (3), since all other variables are known from the 
experimental data. 

The emissivity and absorptivity gas depends the 
temperature the gas and moreover the absorptivity depends 
the temperature the radiating source. has shown 
that both the absorptivity and emissivity gas can approxi- 
mated 


empirical constant (equal 0.45) used evaluating gas absorp- 
from gas emissivity data. found from the available 
emissivity the average the surface and gas tempera- 
tures (this procedure not recommended when the surface and 
gas temperature differ more than factor two). Using 
these average values gas emissivity, the absorption 
and the fraction may calculated for this case from 
(6) and (7). 

Cohen(7) has shown that the emissivity gas may also 
expressed volume basis and that equal 4k. Hence 
the emission rate gas volume 


E= 4kdV o 


Returning the main problem, the reception factors between 
any two zones the model must known function the 
gas absorption strength before solution can attempted. 
general, these factors express the fraction the energy originat- 
ing any zone which reaches and absorbed 


more accurate estimate the mean beam length 
enclosure is required, use may be made of the following relationship 
(based on an analysis presented in a recent paper(!”)): 
where 1.0, 1.25, 1.5, 2.0, 2.5, 3.0, 4.0. 
and the mean beam length, defined the equation: 


4k V, ®,/A, =1-—e 
The correction in the present case is negligible. 
(8) has indicated that the emission rate per unit area 


of black surface into a vacuum and that this should be modified for 


other media by multiplying by the square of the refractive index. Since 
the refractive index of most gases is near one, this factor can be ignored. 


Figure 2—Reception factors for spher- 
ical system Figure for 
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Figure 3—Reception factors for spher- 
ical system Figure for 


another zone. Surfaces are considered completely absorb- 
ing. The formulation these factors given Appendix 
and the resulting these factors function are 


has been termed direct-interchange and 


designated Similarly definition 


and 


since, from the first two laws thermodynamics, 
interchange must zero when the temperatures these zones 
are the same. For gas-to-gas interchange 


When using these interchange areas, each must multiplied 
the fraction since the gas capable radiating and absorbing 
energy only over fraction the wave- Over the re- 
mainder the spectrum the gas-to-gas and gas-to-surface factors 
are zero. 


These factors are applicable black-walled system; 
those cases where the droplet does not have high absorptivity, 
some consideration must given the transmitted and reflected 
part the radiation striking the drop. These effects can 
accounted for approximately the methods described 
Appendix 


Since sharp temperature variations occur the boundary 
layer, effective temperature must assigned each gas 
zone. The effective temperature any gas zone must account 
for the increasing volume gas points furthest removed 
from the drop, well for the changing temperature over the 
given zone. definition, the effective temperature must 
such that: 


= 4ka Vi; o 
ri 


The subscripts and indicate the inner and outer limits the 
gas zone. Over the small radial distances question, may 
considered vary linearly with such that: 
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This allows Equation (15) integrated yield: Expressed 
t t 
t 


where 


The ‘al bal Equation (3) can used calculate approximate 
and balances each gas zone, thus temperature profile the boundary layer, which provides 
take the convenient starting point. this way, temperatures may 
assigned the boundaries each zone and the effective 
Heat conducted into Heat conducted out temperatures calculated from Equation (17). The conduction 

zone from the outer zone into the adja- (2) 
adjacent zone. appear unknowns the balance equation. However, 
the evaporation rate the drop surface (subtracting the direct 
Sensible heat associated Sensible radiation contribution) boundary condition which enables 
with vapor flow into with the vapor flow out these terms evaluated every point the boundary 
from the inner adjacent gj. 
layer (Equation (19)). Substituting the values the thermal 
conductivity and area the respective boundary allows 
Radiation into from Radiation leaving calculation the temperature gradient each point. The 
all other surface and gas all other surface and gas resulting curve (dT/dr),, versus can graphically inte- 


grated yield the temperature distribution which accounts for 
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Figure 7—Reception factors for spher- Figure 8—Reception factors for spher- Figure 9—Reception factors for spher- 
ical system Figure for ical system Figure for ical system Figure for s,. 
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TABLE 


RESULTS CALCULATED TEMPERATURE PROFILE BOUNDARY 
WITH AND WITHOUT RADIATION THE GAS 


10-4 (inner 
Gas Zone (No Radia- (Radiation 
(With and 
Without Radiation) 
1.00 48.8 48.5 1.497 100 
29.9 0.772 235 
2.00 26.1 25.4 0.256 402 
2.50 10.6 11.9 0.158 453 
3.00 0.110 487 


the radiation absorption phenomenon and changing gas properties 
the boundary layer. This iterative procedure repeated until 
the calculated temperature profile coincides with the assumed 
one. those cases where the gas absorbs appreciable amount 
radiation, the direct radiation the drop will have 
reduced accordingly, which turn results higher tempera- 
ture gradient the drop surface effect the same evaporation. 
This calculation can also performed excluding the adsorp- 
tion/emission radiation the gas, but including the effect 
temperature thermal conductivity and heat capacity. The 
difference then measure the effect radiation 


must pointed out that this calculation not, intended 
show the effect radiation the boundary conditions, 
evaporation rate and boundary layer thickness. merely 
gives appreciation the magnitude this phenomenon and 
its effect the temperature profile under given set condi- 
tions and for the assumed model: 


Example 


Diameter water droplet 1000 microns 
Temperature gas and wall 500°C. 

Evaporation rate 73.5 
Direct radiation the droplet 

without attenuation the 

boundary layer (a, 0.90) 12.9 
[Equation (3)| 


The boundary layer was assumed extend three times the 
droplet diameter for the calculation the mean beam length for 


WATER DROP ~d=10 mm. 
© CALCULATED POINTS 


(RADIATION GAS 
INCLUDED) 
CALCULATED CURVE 
(RADIATION 
NEGLECTED) 


CALCULATED EQN.(3) 


TEMPERATURE (°C) 


Figure 10—Comparison calculated 
results. 


Figure 11— Spherical system upon 
which the formulation f,,,, based. 


the system. the methods already described, and were 
calculated 0.065 and 0.10 respectively. The initial assumed 
temperature profile, shown Figure 10, was calculated with 
and 0.480 (cal.)/(g.) (°C.). The results the calcula- 
tions are shown Table and Figure illustrates the compar- 
ison the temperature profiles with and without radiation 
the gas. 


Discussion Results 


The absorption radiation the boundary layer very 
small and has negligible effect the 
around the drop, and hence does not influence the droplet 
evaporation rate. other words, the gas not sufficiently 
opaque atmospheric pressure absorb appreciable amount 
radiation. Increasing the pressure the gas known 
broaden the absorption and general increase the 
absorption power gas considerably. Under these conditions 
higher pressures and possibly higher temperatures, the 
radiation energy absorbed the boundary layer might become 
which case the temperature gradients around 
the drop would different. This would fact result, since 
the radiant heat transfer direct the drop would reduced, 
and for any given evaporation rate this would necessitate 
sharper temperature gradient the drop surface. points 
removed from the drop surface, the gradients would smaller 
than expected because part the heat load would supplied 
radiation, Equation (19) shows. 


Conclusions 


analysis based spherical symmetrical boundary 
layer model and neglecting any convection effects 
presented for evaporating droplet high temperature sur- 
roundings. iterative method solution was adopted, based 
mass and energy balances which took into account heat 
transfer conduction and radiation well the sensible heat 
transfer associated with radial vapor flow, under conditions 
changing gas properties with temperature. was considered 
real, but its absorption strength was assumed constant 
throughout the spherical shell. should noted that this 
method solution quite general and its application 
cylindrical system has been discussed recent 

Calculations showed that for drop about one millimeter 
diameter, the temperature profiles with and without absorption 
radiation the boundary layer are practically identical. 
These findings substantiate the assumption negligible absorp- 
tion radiation the boundary layer which was made 
previous impossible, the present state 
extend this analysis take into account the 
convection currents set the cold vapor the boundary 
layer evaporating droplet. Although these convective 


effects will undoubtedly modify the temperature profiles calcu- 
lated the present analysis, they could not conceivably cause 
greater absorption radiation the boundary layer. 


Figure Spherical system upon 
which the formulation based. 
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Figure 13—Spherical system upon which the formulation 
based. 


important aspect this study was the calculations 
accurate reception factors for spherical systems radius 
ratio (r/ra) which are presented Figures Total 
interchange areas have also been worked out (see Appendix 
which allow for reflection the surface the inner sphere. 
believed that this information will become increasingly 
useful, view the growing number high- 
and applications which radiant heat transmission 
plays important part. With slight modifications, the analysis 
and method solution presented this paper could applied 
the evaporation and combustion sprays liquid fuels, the 
reduction particles metallic oxides, and various other 
chemical reactions particulate systems. Although based 
model, the information thus obtained would provide 
better understanding the actual, more complex situation. 

interesting note that very recently analysis 
involving heat-generating gas spherical system similar 
the one discussed this paper was presented Sparrow, 
Usiskin and neglecting conduction and convec- 
tion, and assuming that the internal rate heat generation 
was uniform, that the gas was grey and uniform absorption 
coefficient (independent temperature), and that the spherical 
enclosure was black-walled, energy balance infinitesimal 
spherical volume was established, which permitted the in- 
genious formulation governing integral The latter 
was solved numerical means using iterative method and 
the results for enclosing surfaces equal temperature were 
given the form plots the gas emissive power distribution 
(or temperature distribution) the enclosure, function 
the geometry size the system and for wide variation 
the value the absorption coefficient. believed that the 
general method solution discussed the present paper could 
have been applied this special case; its use would probably 
required for the solution the more general case, where the 
two enclosing surfaces are different temperatures. 
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Nomenclature 


Roman Symbols 

weighting factor Equation (5) 
surface area, cm.? 
total surface boundaries enclosing volume gas 
exponent temperature ratio which allows the gas 
absorptivity evaluated from emissivity data 
heat capacity, cal./(g.) (°C.) 
diameter effective spherical boundary layer, cm. 
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Figure 14—Spherical system upon which the formulation 
based. 


radiant energy emission rate, 

direct geometric view factor for black-walled system 
containing non-absorbing medium fraction radia- 
tion leaving one zone which reaches and absorbed 
another (subscripts denote sender and receiver) 


reception factor for black-walled enclosure contain- 


ing absorbing medium fraction radiation leaving 
one zone which reaches and absorbed another 
(subscripts denote sender, first, and receiver, second). 


fraction radiation leaving elemental volume gas 


which reaches and absorbed inner spherical 
black surface 
fraction radiation leaving elemental volume gas 
which reaches and absorbed outer spherical 
black surface 
gas zones 
total interchange areas for grey-walled system 
containing absorbing medium area regulates the 
amount radiation leaving any zone (first letter) which 
reaches and absorbed another zone (second letter), 
including the reflection surface zones, cm.? 
direct interchange areas for black-walled system con- 
taining absorbing medium, calculated from f-factors, 
cm.? 
subscripts denoting any zone 
absorption coefficient for the gas, cm.~! 
mean beam length, cm. 
thickness absorbing medium, cm. 
subscripts denoting any zone 
mass evaporated, 
partial pressure absorbing gas, atm. 
radiant heat transfer rate droplet, cal./sec. 
radiant heat transfer rate gas between droplet surface 
and radius cal./(sec.) (unit volume) 
separating distance, cm., unless normalized 
radius any point effective boundary layer, cm., 
unless normalized 
radius effective spherical boundary layer, cm. 
droplet radius, cm. 
surface zones 
absolute temperature, °K. 
droplet temperature, °K. 
gas temperature, °K. 
gas volume, total gas volume; volume 
defined Equation (20) 


Greek Symbols 
absorptivity 
absorptivity particle 
absorptivity surface 
constant thermal conductivity equation 
emissivity 
emissivity gas 
emissivity evaluated one and two mean beam lengths, 
respectively 
surface 
defined Equation (18) 
time, sec. 
angles Figure 11, 
latent heat vaporization, cal./g. 
absolute viscosity, g./(sec.) (cm.) 
total pressure, atm. 


ith 
la- 
/ | 

ile 
let 
ly 
[- 
r 
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reflectivity 
Stefan-Boltzmann constant: 1.355 (cal.)/(sec.) 
(cm.?) 
spherical co-ordinate, Equation (1); also angles 
Figures and 
escape factor 
spherical co-ordinate, Equation (1) 
factor which characteristic length multiplied 
give the mean beam-length the system 

solid angle, Figures 11, 12, and 
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APPENDIX 


Formulation the Reception Factors for 
Spherical System 


the analysis the radiation absorbed the boundary 
layer evaporating drop, spherical symmetrical model was 
postulated. The region around the drop which the temperature 
gradients exist was divided into number spherical gas shells, 
that the distribution absorbed radiation each gas shell 
could calculated (see Figure 1). 

formulate the reception factors for this system, two 
assumptions must made: 

(1) the surfaces are black 

(2) the absorption strength the gas uniform throughout 

The latter assumption has been discussed detail Hottel 
and 

The following types reception factors are considered: 

(1) surface-to-surface 

(2) gas-to-inner and -outer surfaces 

(3) gas-to-gas 

make the solution completely general, all length terms are 
normalized with respect the droplet radius (inner sphere) 


(1) Surface-to-surface. 
(a) Inner outer surface. 


Referring Figure 11, consider the radiation leaving 
angle the normal The quantity radiation received 
dA, located distance proportional the relative area 
each seen the other and inversely proportional the 
square the separating distance. Moreover this radiation will 
reduced result absorption the intervening gas. 
Introducing the proportionality constant the one way 
radiation from becomes: 


definition 

Equating and (A2), and since dw, 


TAs, A, 


Since every part A,, sees the same relative portion 


must independent where located and hence: 


= am | da 


Equation (A4) becomes: 


where the limits integration are 


Then 


where 


Both integrals can evaluated integrating parts, hence: 


b 
— 


ponential integral, which has been and and 
are given (A6). 


(b) Outer surface outer surface. 
Referring Figure 12, the relation which describes the fraction 


the radiation leaving the solid angle and intercepted 
the outer sphere given by: 


which can written (since 


Integrating parts: 


0 


(2) Gas-to-surface factors. 


(a) Gas the inner surface sphere surface. 


Referring Figure 13, volume element dV, located point 
radiating over the total solid angle steradians. The radiation 
reaching dA, located distance from contained the 
differential solid angle dw. This radiation must reduced 
the transmissivity the gas over the separating distance. Hence, 
the one way flux expressed by: 


—kraR 


Substituting for sin® from Figure 13, dividing the 
total emission from volume element and integrating over the 
solid angle which includes the inner surface, the result the 


the system considered as 
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absorbed 
as at p 


Equation 


kra[- 
where 


(b) 
Referring 
from 
which 


Dividing 
grating 
yields 
originatin 


Substitut 


where 
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The 
(whether 


Since 


The num 
methods. 


(3) Gas 
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from the 
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the radiation originating which reaches and 
absorbed the inner spherical surface. This fraction defined 


1 


Equation (A15) can integrated parts, twice, yield: 


(A17) 


(b) Gas any outer surface. 


Referring Figure 14, the equation describing the radiation 

from volume element contained differential solid angle 

which reaches and absorbed the outer sphere located 
is: 


Dividing the total emission volume element and inte- 
over that portion the outer surface which ‘‘sees”, 

yields the following expression for the fraction the energy 

originating which absorbed the outer surface: 


Substituting for sin® from Figure 14, and integrating: 


(1/kra) fe 


The gas-to-surface factors can found from the 


Vv 


Since the same for all points given radius 


Va, 


The numerator (A23) conveniently integrated numerical 


methods. 


(3) Gas-to-gas factors. 


The reception factors for any two gas zones may obtained 
from the expressions. The difference between the factors for 
any two imaginary surfaces bounding the gas zone will yield the 
fraction the radiation originating which absorbed 
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the given volume. Integrating over the emitting zone volume 
and dividing the total emission from that gas zone yields the 
desired gas-to-gas reception factor, 

Thus 


Ve; 


The integral can evaluated numerical methods. 


Evaluation the Factors 


The expressions for the factors were evaluated means 
650 digital computer. The numerical integrations 
were carried out the computer, employing Simpson’s Rule. 

The integrals were evaluated desk calculator. all 
cases, eight place accuracy was maintained throughout. 

The final factors are indicated Figures and Table 2.* 


APPENDIX 
Allowance for Reflection the Drop Surface 


The direct interchange areas are obtained from the reception 
factors. However, these areas are applicable only black 
system. imaginary black surface can postulated the 
outer edge the boundary layer; the inner drop, however, may 
only partially absorbing and some the radiation striking 
this surface may reflected and some transmitted. The problem 
becomes even more complicated attempt made trace 
the transmitted beam. However, one assumes the transmis- 
will small, since the radiation which not absorbed the 
drop still being considered. If, moreover, the surface assumed 
act grey body, the absorption the reflected beam can 
easily calculated from the direct interchange areas. account 
for this reflected flux, another interchange area introduced, 
which corresponds (AF) the old These 
terms are known total interchange areas and are defined 
follows: 


GiG; the fraction the energy originating (si, 
respectively) which absorbed respectively) taking 
into account the direct radiation, well the multiple reflections 
all surfaces the enclosure, and the absorption the inter- 
vening gas. 

The determinant method for formulating these total inter- 
change areas from the direct interchange areas was discussed 
recent paper(*). Applying this method the present case, the 
following total interchange areas result: 


(1) 


When the absorptivity the droplet high (a,, 0.90), the 


increased absorption the gas, arising result reflection 
the drop surface, small and may neglected without 
seriously affecting the result. 


this paper has been deposited Document No. 6878 with 
the ADI Auxiliary Publications Project, Photoduplication Service, Library 
of Congress, Washington 25, D.C. A copy may be obtained by citing 
the Document No. and by remitting $1.25 for photoprints, or $1.25 for 
35 mm. microfilm. Advance payment is required. Make cheques or 
money orders payable to: Chief, Photoduplication Service, Library of 


Congress. 
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Sulphur and Sulphur Dioxide for 
Separating Copper from 
Nickel-Cobalt-Copper Sulphate 


hydro-metallurgical processes for the produc- 
tion nickel and cobalt, the leach solutions usually 
contains salts nickel, cobalt and copper. recover 
high purity nickel and cobalt metals from such solu- 
tions, the copper removed prior the nickel-cobalt 
separation. 

During the course extensive 
determine the most economical and efficient methods 
for separating the copper, was established that 
copper solution reacts with elemental sulphur 
the presence sulphur dioxide forming copper 
sulphide precipitate. This reaction permits simple 
and quantitative separation copper 
and cobalt leach solutions. 

The variables affecting this reaction were studied 
and the results are presented this paper. Studies 
were conducted on: the sulphur concentration be- 
tween 2.0/1 molar ratio sulphur copper, the 
sulphur dioxide concentration gpl., the tem- 
perature the range from 100°F. 325°F., the 
copper concentration from gpl. and from 
1.0 9.0. 

The results are correlated and tentative mech- 
anism for the reaction discussed. 


and cobalt concentrates mattes prepared from 
nickel-cobalt-sulphide ores usually contain 
quantity copper. Among number conventional methods 
for processing these oxidizing leaches either 
ammoniacal ammonium sulphate solutions sulphuric acid. 

During the above leach treatment, the nickel, cobalt and 
copper sulphides are converted water soluble sulphates and 
are separated from the leach residue 

The compositions the leach solutions, each case are 
dependent upon the composition the treated material. 


The compositions typical leach solutions are shown below: 


Type Ni | Co | Cu | NH3F | (NH4)2S0j | H2SO4 
of Solution | gpl. | gpl. |} gpl. | gpl. gpl. | gpl. pH 
Ammoniacal | 50 0.5-2.0 | 5-10 | 30-100 | 100-300 0 | 7-9 


\Manuscript received March 27; accepted July 27, 1961. 

Research and Development Division, Sherritt Gordon Mines Limited, 
Fort Saskatchewan, Alta. 

Based on a paper presented to the Joint A.I.Ch.E.-C.I.C. Chemical Engin- 
eering Conference, Cleveland, Ohio, May 7-10, 1961. 
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Copper must removed from these solutions before the 
recovery nickel and cobalt reduction with hydrogen 
clevated temperature and 


has been discovered that this removal copper 
achieved the addition sulphur the leach 
solutions thereby causing the precipitation copper sulphide. 
This reaction proceeds only the presence sulphite, thio- 
sulphate and polythionate tons. 

The ammoniacal leach solutions may contain 
amount these compounds for the precipitation copper 
sulphide accomplished the addition sulphur alone. 
However, the concentrations thionates acid leach solutions 
are too small for reasonable rate precipitation and 
therefore advantageous add sulphur dioxide source 
sulphite, thiosulphate and thionate 

This paper describes the conditions under which copper can 
precipitated from sulphate solutions copper sulphide 
the addition sulphur and sulphur dioxide. 


Experimental Procedure and Results 

The solutions used this work were either actual acid 
ammoniacal leach solutions obtained 
synthetic solutions with the same general composition leach 
solutions. 

The experiments were carried out 
steel autoclave equipped with mechanical agitator and heated 
cooled externally gas water respectiv 

The autoclave was charged with litres solution and 
calculated quantity sulphur powder. dioxide 
ammonium sulphite solution was added usually the reaction 
temperature either directly through the charging port 
means pressure bomb the temperature was ‘above 

was determined the head solution titration with 
0.008 solution. The same titration was used during 
the reaction determine combined sulphide, thiosulphate and 
polythionate concentrations. However certain specific tests 
the sulphite, thiosulphate, polythionate and sulphate concentra- 
tion were determined separately (see Figure 14). The sulphite 
concentrations were determined from the differences titration 
equal volumes the sample acetic acid and formaldehyde 
solutions using standard todine solution. The thiosulphate 
concentration was determined from direct 
iodine the formaldehyde solution. The polythionates were 
determined from the differences between the iodine titration 
the acetic acid solution and third sample hydrochloric acid 
solution using standard potassium bromate solution. 
was used the indicator for the iodine titrations and methyl- 
orange for the titration. 
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Figure 1—Effect varying the sul- 

phur addition the rate copper 

precipitation from ammoniacal Ni- 
Co-Cu sulphate solution. 
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Figure 2—Effect varying the sul- 

phur addition the rate copper 

precipitation from acid Ni-Co-Cu 
sulphate solution. 


10.0 


T= 
SO,= O GPL. 
x x x x x x 
> 
z 2.2GPL 
z 
22° 7.6 GPL 
= 
z 
° 
9.9 
wi 195° F 
so, = VARIABLE 
cu = 4,6 GPL 


10 20 30 40 SO 60 70 80 90 


0,001 


RETENTION TIME (MIN) 


Figure 3—Effect sulphur dioxide 

concentration the rate copper 

precipitation from ammoniacal Ni- 
Co-Cu sulphate solution. 


The sulphate ion concentration was determined gravimetri- 
cally direct precipitation barium sulphate using barium 
chloride. The copper concentration solution was followed 
taking samples suitable intervals and analyzing for copper 
electroplating colorimetrically. 

After completion the test the precipitated copper sulphide 
was recovered, washed and analyzed. 

During the preliminary experiments was observed that the 
extent and rate copper sulphide precipitation was influenced 
the concentrations sulphur, sulphur dioxide and copper 
and also the temperature and the hydrogen ion concentration 
the reaction system. 

The concentration nickel was found without effect 
the range 60-120 gpl. acid solutions and 40-50 gpl. 
ammoniacal solutions. The same holds true for the ammonium 
sulphate concentration 100 gpl. acid, and between 200 
and 300 gpl. ammoniacal, leach solutions. 


The Effect Sulphur Concentration 


The reaction between sulphur and copper sulphate solu- 
tions heterogeneous reaction and therefore the physical 
characteristics the solid phase are very important. Decreasing 
particle size sulphur and treatment the sulphur with wetting 
agents accelerate the reaction rate. 

Technical reagent grade sulphur passing 200 mesh and 
treated with Aerosol was found very effective for 
the precipitation copper. 

The quantity sulphur added the solution was expressed 
the sulphur copper molar ratio. The relationship between 
the S°/Cu molar ratio and the rate copper precipitation under 
the conditions specified Figures and showed that the rate 
copper precipitation increased, and copper content the 
purified solution decreased, with increase sulphur concen- 
tration. The minimum quantity sulphur required for copper 
precipitation between 0.5 1.0 mols sulphur per mol 
copper solution. 


Effect Sulphur Dioxide Concentration 


The reaction between elemental sulphur and copper ions 
sulphate solutions leading the precipitation copper sulphide 
are complex and involve sulphite, thiosulphate 
species, which may act carriers and reductants the sulphur. 

Leach solutions obtained the extraction sulphide 
highly ammoniacal solutions ammonium 
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sulphate usually contain sufficient concentration these inter- 
mediate species give acceptable rate copper sulphide 
precipitation when elemental sulphur added the leach 
solutions. However, the concentration these sulphur carriers 
acid leach solutions, and some ammonium sulphate leach 
solutions containing relatively small amounts ammonia, too 
low give acceptable rates copper sulphide precipitation 
when sulphur powder added these solutions. these cases 
external source sulphite, thiosulphate thionate must 
added with the elemental sulphur. This addition most easily 
accomplished the addition sulphur dioxide, ammonium 
sulphite. 

The effect the amount sulphur dioxide added shown 
Figures and from which can seen that the rate 
copper sulphide precipitation function the sulphur dioxide 
concentration. 


The Effect Temperature 


The effect temperature was investigated between 100° 
and Figures and show that the temperature has very 
pronounced effect upon the rate copper sulphide precipitation. 


Because this marked temperature-dependence the quantity 
sulphur dioxide required achieve certain rate copper 
sulphide precipitation could decreased increasing the 
reaction temperature. This shown Figures and The 
lowest temperature for satisfactory rate copper sulphide 
precipitation appears 125°F. while above 325°F. copper 
sulphide precipitation occurs even when sulphur dioxide not 
added the solution. 


The Effect Copper Concentration 


The copper concentration leach solutions may vary widely 
depending upon the composition the treated ‘material, but 
usually the copper concentration solution between and 
gpl. These investigations however, were extended beyond 
this limit and the copper precipitation was carried out solu- 
tions containing gpl. copper. The results illustrated 
Figures and show that below 212°F. the initial copper 
concentration has negative effect upon the copper sulphide 
precipitation and more sulphur dioxide required precipitate 
the copper the specification level. The relationship between 
initial copper concentration and sulphur dioxide concentration 
212°F. shown Figures and 11. 

higher temperatures the effect the sulphur dioxide 
concentration less pronounced and 270°F. (Figure 
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ceases completely. This can explained the fact that 


low temperatures the reaction affected the presence 


thionates. Increasing the copper concentration the solution, 
then causes the thionates used more rapidly than they are 
produced and the rate copper precipitation decreases. 
higher temperatures however, the direct reactivity the ele- 
mental sulphur increases and the reaction proceeds completion 
soon the required quantity elemental sulphur available. 


The Effect Hydrogen Ion Concentration 


The effect was studied both the ammoniacal and 
the acid leach solutions. 

Figure illustrates the effect increasing amounts 
ammonia the precipitation copper sulphide from 
ammoniacal leach solution. The amount ammonia present 
expressed molar ratio which the ratio the number 
mols ammonia which can titrated with sulphuric acid 
the number mols metal (Ni, and Cu) present solu- 
tion. can seen that increase this 
molar ratio increases the final equilibrium concentration 
copper and prevents complete precipitation the copper 
copper sulphide. This another example complex-formation 
exhibiting pronounced effect metal ion reaction aqueous 
solution. 

Figure shows the effect ammonia upon the precipitation 
copper sulphide from the acid leach solutions. this case 
the amount ammonia added expressed NH3/Cu 
molar ratio. The theoretical quantity required determined 
from the following equation: 


From this equation can seen that two mols ammonia 
must added for each mol consumed and two mols 
ammonia for each mol precipitated copper sulphide 
order neutralize the acid produced these reactions. The 
SO, molar ratios Figure calculated according 
Equation (1) using solution with initial 1.7 are 
good agreement with this conclusion. The concentration 
used calculating the SO: molar ratio was the 
amount consumed the reaction and not the 
concentration present the initial solution. The addition 
less than the theoretical amount ammonia caused drop 
which prevented the complete precipitation copper 
sulphide. The use the theoretical amount ammonia 
allows the precipitation proceed completion while the 
remains essentially unchanged. excess ammonia added 
the rises rapidly and above about 4.5 the copper precipi- 
tates out the hydroxide rather than the sulphide. 
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Figure 14—Changes concentration copper and sul- 

phur-dioxide-thiosulphates-trithionates and sulphate ions 

during precipitation copper from ammoniacal Ni-Co-Cu 
sulphate solution with sulphur and sulphur dioxide. 


THE REACTION PRODUCTS AND REACTION 
MECHANISM 


The Nature the Precipitated Copper Sulphide 

The material precipitated during the reaction high grade 
copper sulphide contaminated with small quantity elemental 
sulphur and nickel sulphide. 

Typical analyses the residues precipitated from the two 
types solution are: 


Cu | Total S |Elemental S| Ni Co Sr-—S°/Cu 


Ammoniacal d 3 0.8 <0.01 0.8 
Acid 0.1 <0.01 1.0 


The copper sulphide precipitated from ammoniacal solution 
characterized higher contamination with nickel sulphide. 
This explained the fact that sulphide ions will precipitate 
both metals above However, the copper preferentially 
precipitated and only small amount nickel sulphide present 
it. 

acid solution sulphide ions not precipitate nickel and 
the copper sulphide therefore contains very little nickel when 
precipitated under these conditions. 

The sulphide sulphur, which was calculated difference 
between total sulphur and elemental sulphur, copper ratio 
the copper sulphide indicates that ammoniacal solutions 
copper precipitates mixture CuS and whereas 
acid solution the precipitate cupric sulphide. 

The elemental sulphur the copper sulphide may deviate 
widely from the typical analysis depending upon the precipitation 
conditions. 


The Reaction Mechanism 

The mechanism the reaction between copper ions and 
elemental sulphur produce copper sulphide complex one 
and involves many competing, simultaneous and consecutive 
reactions, all with their own rate temperature-dependence. 

The postulated sequence reactions follows: 
(a) Dissolution sulphur: 

SO, + H.0 SO;-- + 2H+ 

Sk 
(b) Formation thionic acids: 

2S,0;--+ 8Cut++4H,0 + 8Cu+ + SO,-- + 8H*+ .(5) 
(Reaction (5) only possible alkaline solution) 
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(c) Precipitation copper sulphide: 


The order and rate these reactions depends the pH, 
temperature and concentrations the various intermediate ionic 
species. 

The net result these reactions that low temperature 
elemental sulphur reduced sulphide while oxidized 
sulphate, and high temperature the elemental sulphur 
both reduced sulphide and oxidized sulphite. 


The sulphide ion produced: gradually and reacts with the 
copper solution produce copper sulphide. Since the solu- 
bility products NiS, CuS and are 8.5 
and respectively follows from the lack nickel 
sulphide the precipitate that the concentration sulphide 
ions solution does not normally exceed 
under the reaction conditions. 

During the reaction shown Figure the concentrations 
sulphate and thiosulphate rose while the concentration 
sulphur dioxide fell. The decomposition trithionate ions was 
rapid that could detected only the initial stages the 
both these species could demonstrated 
similar solutions the absence metals. 

the temperature increased all these various reactions 
are accelerated until 325°F. and above the copper and elemental 
sulphur can react directly without the intermediary thionate 
ions and the following reaction appears take place. 


(9) 


Conclusion 


The separation copper from solutions containing 
nickel, cobalt and copper sulphates can achieved adding 
elemental sulphur and sulphur dioxide which cause the precipi- 
tation copper sulphide. 

This reaction involves sulphite, thiosulphate 
intermediates and the primary purpose the sulphur 
dioxide addition generate these ionic species. Above 
these intermediaries, and sulphur dioxide, are unnecessary and 
the addition elemental sulphur sufficient bring about the 
precipitation copper sulphide. 

The reaction can carried out over wide temperature 
range, the lowest practical temperature being 125°F. The 
dependence rate upon temperature very great. 
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The reaction can carried out over the range 
although the rate slower these extremes and the small 
range about neutrality cannot used because the precipita- 
tion copper hydroxide. 


The acid produced the reaction can neutralized the 
base most suitable for the subsequent steps the process. 
may include NaOH, KOH, CaO, MgO, 
and also metals such cobalt and nickel 
can used for acid neutralization. 

Under properly controlled conditions the rate copper 
removal very rapid and copper can removed less than 
0.02 gpl. within few minutes. The precipitated copper sulphide 
high grade material containing less than nickel and cobalt. 

This method copper separation very flexible and can 
adapted various hydrometallurgical processes, where the 
main reagents, and sulphur, are readily available. 
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Research the McGill Hypersonic 
Propulsion Research 


supersonic combustion ramjets started 1958. 


Theoretical studies show that the performance good 
Mach 25; they are thus suitable for satellite boosting 


Facilities which have been constructed for experimental 


work are: 


(1) Zirconia pebble bed heat 0.1 air 2000°C. 
100 
(2) Vacuum system 0.1 0.5 connected 


test facilities. 


(3) Hypersonic shock tunnel, Mach upwards. 

Fuel mixing study facility providing two supersonic 
dissimilar gas streams. 

(5) Nozzle recombination test facility using nitrogen dioxide 
diluted. 

(6) Thrust measuring test facility. 

(7) High temperature electric furnace. 

Most the combustion research date has been carried 
out using air heated the zirconia pebble bed, and has been 
directed obtaining information relevant shock free 
supersonic combustion. 

The following experiments have been carried out. 

(1) Spontaneous combustion hydrogen Mach 

This experiment was simply give the lower ignition 
hydrogen flow system. temperature was 
found about 600°C. 

(2) Ignition delay hydrogen/air mixtures high temper- 
atures. 

This data was acquired extend the data currently available 
1300°K. sub-atmospheric pressures. The results are 
shown Figure 

(3) Boundary layer 
and 2.6. 

These experiments were carried out determine the possi- 
bility supersonic combustion with porous wall fuel injection. 
Difficulties were encountered due separation the boundary 
layer caused combustion heat release. This considered 
inherent with this type fuel injection, and has been 
discontinued. 

(+) Downstream injection hydrogen 
slotted nozzle. 


1.2 using 


Satisfactory combustion was obtained with change 
fame character the Mach No. went through The length 
diameter ratio the flame retained its classical relationship 


‘Manuscript received September 7, 1961. 
‘Mechanical Engineering Department, McGill University, Montreal, Que. 
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throughout. Higher Mach Numbers and/or lower heat release 
was required for shock combustion. 

(5) Downstream and cross-stream injection tri-ethyl 

Using air with reduced oxygen content decrease the 
maximum heat release showed that combustion could easily 
obtained, although certain precautions had taken with 
the injector prevent plugging due decomposition the 
T.E.A 

cross-stream injection. 

This experiment used radial stub jets penetrate the 
boundary layer; the combustion was carried out instru- 
mented constant area tube. The results showed that shock 
formed near the fuel injector station and great difficulty was 
caused the aluminum oxide deposit fouling the instrumenta- 
tion lines. 

(7) Constant area combusion hydrogen using 
cross-stream injection. 

These experiments are present underway, using air total 
temperatures 2500°C. obtained vitiating the air. 


Future supersonic combustion research requires even higher 
temperatures. However, since the ignition delays are less than 
seconds, believe that satisfactory results will obtained 
intermittent type facilities with running times few 
milliseconds. The proposals are discussed follows: 


LEGEND 
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Figure 1—Relation versus constant pressure. 
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Studies Hypersonic Air-Breathing Engines 
the Shock Tunnel 
Combustion Research 

high pressure (80,000 p.s.i.) shock tunnel being con- 
structed study combustion phenomena conditions occurring 
supersonic combustion ramjet 


Total 
Total pressure 76,000 
Mach Number 


These conditions will simulated the nozzle exit the 
shock tunnel. The initial experiments will involve downstream 
injection through axial injector through the main shock 
tunnel nozzle. Subsequent tests will show the effectiveness 
cross-stream injection. The fuel, hydrogen most cases, will 
injected into the main airflow means auxiliary shock 
tube, driven either its own driver from the main testing 
air driver. 


Intake Research 


Experimentation and theoretical work progress study 
the effects shock waves and boundary layers the efficient 
deceleration hypersonic flow supersonic speeds. Straight 
wedges and oblique shock models are presently being tested 
12. Future experiments will involve curved intake surface 
with controlled boundary layers. 


266 


Nozzle Research 


Both the high and the low pressure shock tube will used 
for nozzle flow studies. Combustion products the driven 
tube will shock heated/compressed and blown into propul- 
sion nozzle. Pressure and thrust measurements will then deter- 
mine the performance. 


Connected Tests 

After completion the above tests together with the present 
intake tests planned run experiments complete con- 
figurations, i.e., intake, combustion and exhaust. Subsequent 
experiments will then conducted with gun-launched hyper- 
sonic free-flight models. 
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INDUSTRIAL SECTION 


The Design and Operation 
High Pressure Reforming Pilot Plant’ 


pilot plant has been designed and constructed 
for the exploration reforming processes temper- 
atures 1200°F. and pressures 1500 p.s.i.g. 
The components the pilot plant are described, 
well the calculation material balances and pro- 
duct yields IBM 650 computer. 


hen catalytic reforming was introduced Imperial Oil’s 

refineries, was decided build reforming pilot plant. 
The pilot plant was evaluate the effects new feeds, catalysts 
and changes the operating conditions, well assist the 
refineries any problems that may develop. This paper presents 
brief outline the final design, incorporating some modifica- 
tions carried out following the start- operation. 


Design Basis 


The first problem solved was the size the pilot unit. 
For simplicity operation, and minimize catalyst and feed 
requirements, the unit had small possible, consistent 
with obtaining adequate sample for quality evaluation 
reasonable time. But another factor had taken into account. 
Many reformer feeds contain more sulphur than can 
the reforming catalyst, and they have hydrogen treated 
hydrofined before reforming. Normally this would not 
important, since the feed could obtained from the refinery, 
after hydrofining. However, the discovery new oil fields 
presents new naphthas, not processed refinery, for evalua- 
tion. These would have drofined the laboratory before 
Since the drofining flow plan not very different 
from that reforming, the same pilot plant could used, but 
should have relativ ely large capacity, that too much time 
would not spent feed preparation. the final design, this 
problem was solved the use interchangeable reactors 
different capacity. The rest the system was sized process 
the expected maximum feed rate. 


Catalytic reforming carried out 850-950°F. and 300-500 
Over platinum catalyst improve the octane rating 
gasoline conversion octane components. Figure 
shows the principal reactions involved. The most desirable 
teactions are the conversion naphthenes 
aromatics. Both these reactions produce high octane material 
ely low octane components, and are endothermic. 
straight chain slightly exothermic, 
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while hydrocracking fairly exothermic. The over- 
all reaction usually quite endothermic, causing appreciable 
temperature drop the reactor. 


its simplest form, the catalytic reforming pilot plant flow 
plan can represented the scheme shown Figure The 
feed, consisting mainly naphthenes and paraffins taken from 
storage, mixed with hydrogen recycle gas, and passed into the 
reactor. The product enters separator, where hydrogen rich 
gas removed from the hydrocarbons, and then the aromatics 
rich reformate stabilized removing the light hydrocarbons— 
ethane, propane and butanes. The hydrogen rich gas recycled. 
Provision has been made introduce other gases drogen, 
nitrogen air into the system, shown the dotted section 
Figure 


Reactor Design 


The heart the entire system the reactor section. Here 
the incoming feed has heated the reaction temperature, 
and contacted with the catalyst. 


There are two possible ways operating the reactor 
isothermally adiabatically. the isothermal operation enough 
heat supplied the reactor maintain constant temperature 
throughout the catalyst bed. This method would preferable, 
since simplifies the evaluation the reaction mechanisms 
involved. adiabatic operation, only enough heat supplied 
the reactor balance heat losses. Here temperature gradient 
may exist along the reactor. This approximates commercial 
reactor, and thus may more indicative the performance 
expected. Figure shows the large difference the tempera- 
ture profile between adiabatic and isothermal reactor, using 
the same feed, catalyst and inlet temperature. The isothermal 
system was not truly isothermal with 75°F. temperature drop 
the first part the reactor, where most the highly endo- 
thermic naphthene dehydrogenation occurred. The temperature 
the sandbath, used heat the isothermal reactor, shown 
for comparison. 

was decided use both methods reactor heating, and 
two parallel reactor systems were designed. The isothermal 
system shown Figure consists single reactor, 
fluidized sandbath. The outside the sandbath 
heated resistance wire heater, and the sand acts the heat 
transfer medium. The fluidizing air preheated aluminum 
block heater, while the feed preheated the coiled part 
the inlet pipe. Two reactors were built, identical dimensions, 
except for the diameter the catalyst chamber. This permits 
easy change the reactor capacity. The temperature the 
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Figure 1—Typical reforming reactions. 
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Figure reactor system. 


reactor controlled thermocouple immersed the sandbath, 
connected controller set hold given temperature. 

The adiabatic system, shown Figure consists four 
reactors connected series. The feed passes through 
aluminum block preheater, and heated the transfer lines 
between the reactors. The line heaters are made two 
“Calrod” heating placed alongside the connecting pipe, 
with all three encased pipe filled with copper shot for better 
heat transfer. 


maintain adiabatic conditions, the reactors are well 
insulated, and heaters are also provided. These heaters supply 
only enough heat prevent heat loss. Therefore, the tempera- 
ture controllers are set maintain given temperature difference 
between the catalyst bed and the insulation, rather than hold 
given temperature. Usually the insulation held the same 
temperature the catalyst bed, but occasionally may higher 
compensate for heat transfer along the reactor wall. 
Figure showed, there large temperature drop along the 
reactor. single heater were used compensate the heat 
losses, too much heat would supplied the bottom part, 
and not enough the top. Therefore, several heaters were 
provided, three for the first reactor, with the largest temperature 
drop, and two for each the others, each with separate 
differential temperature controller. 


Figure shows the method insulating the adiabatic re- 
actors. There are three layers insulation. The inner layer 
loose insulation, contained sheet metal jacket with 
removable bottom cover plate. This insulation can dumped, 
remove the reactor. The middle layer insulation covered 
another sheet metal jacket, which supports the wire heaters, 
and thermowells for measurement the insulation temperature. 
The heaters are covered with layer asbestos tape, and 
layer high temperature insulation. 

change catalyst the reactor, the insulating cap can 
removed, and catalyst discharged opening the end union. 
similar closure the top permits loading catalyst without 
removal the reactor from the system. The same closure 
also used the isothermal reactors. permits rapid change 
the catalyst charge. 

The reactors operate rather severe conditions 
and 1500 p.s.i.g. atmosphere consisting mainly 
hydrogen. Type 347 stainless steel stabilized with columbium 
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Figure 5—Arrangement the adiabatic system. 
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Figure 6—Insulation adiabatic reactors. 


was chosen for the reactors, avoid carbide precipitation, 
particularly when back welding fittings. Subsequent experience 
other units has shown that carbide precipitation apparently 
not serious, and type 316 stainless steel suitable for this 
service. Another problem found with these conditions 
find valve that will stay pressure tight for any length time 
high temperatures. result, the reactor outlet block 
valve had placed after the product cooler. 


these two systems, the isothermal simpler operate, 
and has been used more frequently than the adiabatic system, 
since equivalent temperature the nature the reactor seems 
have little effect. The maximum catalyst capacity 2000 
for each system, and smaller loads can obtained dilution 
the catalyst with inert material. 


Feed and Product Recovery 


The feed system shown Figure consists three 
parts container for the feed, means measuring the feed 
rate, and way getting the feed into the reactor. The critical 
point this system accurate measurement the feed rate. 
There are two ways measuring the feed rate weight 
volume. finally adopted, the volume feed measured 
determining the change the level calibrated blowcase. 
Two pairs are used, and liter capacity, 
because the range required feed rates. diaphragm 
pumps driven common motor were chosen. One pump 
normally used for feed rates about 1500 cc./hr., and the 
other used higher feed rates, 8000 cc./hr. While the 
pumps are the constant displacement type, the flow has been 
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Figure 7—Feed system. 


too variable for their use directly calibrated metering pumps. 
All the feed tanks are nitrogen blanketed and nitrogen pressure 
used transfer the feed from the drum the blowcases 

gas feed used, for example once-through 
itis taken from compressed cylinder, measured rotameter, 
and passed into the reactor. recycle gas used, taken 
from the separator, the flow rate measured rotameter, 
compressed and returned the reactor. handle the 
wide range recycle gas flow rates, diaphragm compressor 
with three pumping heads driven common motor was 
chosen. The heads are connected parallel, and can bypassed 
vary the capacity. Usually two heads are used, with the third 
passed, and used one the others fails. 


Figure shows the make gas system, including the gas 
sampling. The reaction products pass through the cooler into 
pressure separator. The liquid from the separator 
withdrawn through pneumatic control valve the debutanizer. 
Some the drogen rich gas taken for recycle, and the rest 
passed through two back pressure regulators series. The 
first regulator low ers the gas pressure from operating pressure 
100 this point, the gas enters timer-controlled 
solenoid valve, which can pass the gas into sample bombs, 
wet test meter. these lines has another back 
pressure regulator, reducing the pressure sample bomb pres- 
sure the one case, and atmospheric pressure the other. 
This arrangement permits maximum sample bomb pressure 
100 p.s.i.g. The gas going the wet test meter passes through 
prevent evaporation water the meter. 


The debutanizer system similar the make gas system, 
but without high-pressure back pressure regulator. The 
consists electrically heated bottom section, 
and steel pipe filled with stainless steel packing for the 
distillation column. bottoms product withdrawn through 
acontrol valve and cooler into tared sample can. 


Instrumentation and General 

Type 316 stainless steel was used for all vessels operating 
unit pressure, and stainless steel tubing was used for all the 
lines carrying feed and products, regardless pressure. Stainless 
steel was used provide high temperature strength and corrosion 
avoiding contamination products. Carbon steel 
was used for low pressure service, such feed and 
the stabilizer. 

The capacity the pilot plant has presented 
the selection control instruments. There serious problem 
with the temperature control. Conventional controllers can 
used, with thermocouples the sensing elements. liquid 
levels the separator and stabilizer are maintained pneumatic 
control valves, using dielectric probes the sensing elements, 
and pressures back- -pressure regulators. However, the flows 
have controlled manually adjustment feed pump for 
liquid flow and needle valves rotameter inlets for gas flows. 

Since the materials processed hydrogen 
carbons can form explosive mixtures very readily, the entire 
pilot plant had explosion-proof construction. Explosion- 
proof motors, solenoids and conduits for the electrical wiring 
were used, but the use explosion-proof housings 
temperature controllers, electrical switches, timers and auto- 


The Canadian Journal Chemical Engineering, December, 1961 


wer 
SATURAT 
LIQUID 


CONTROL VALVE 8.P.R.-I: BACK PRESSURE 


RE! SET AT 
SURE 
B.P.R.-2: Back PRESSURE 
REGULATOR, SET AT 
100 PSIG 


Figure 8—Make gas system. 


transformers would have been very costly. Therefore, all these 
instruments were housed two air-purged, fully enclosed panel 
boards. boards are maintained under positive air 
pressure all times while the electrical circuits are energized. 
explosivity meter used monitor the unit selected 
points, and give warning flammable vapors escape the 
atmosphere. 


All the high pressure sections are protected safety valves 
connected the vent, and several vent valves are provided for 
the other sections. All pressure vessels have been hydraulically 
pressure-tested 8300 p.s.i.g. shut down the unit, the 
case emergency, such fire, controls have been installed 
outside the building housing the pilot plant. These consist 
two switches cut off the electricity, and blow-down tank 
and valves dump the contents the unit and purge with 
nitrogen. While this installation has not been used for actual 
emergency, used intervals shut down the pilot plant 
the end run, test the equipment and familiarize the 
operators. 

When the actual assembly the equipment was considered, 
there were two choices. One was group all the vessels 
closely possible, minimize the connecting tubing, and thus 
the hold-up the unit. The other alternative, and the one 
actually chosen was place all the vessels line, with 
casy access from both sides. The equipment supported 

“Unistrut” framing two parallel structures. The front 
the unit contains the feed and product recovery sections, all 
the pressure gauges, regulators and most the controlling 
valves. The back section contains the reactors. While this 
assembly involves some long runs piping, has been com- 
pensated the ease with which maintenance and addition 
new equipment can carried out. 

Figure shows overall view the pilot plant. top 
the isothermal reactor and sandbath can seen the centre 
the right-hand window. The adiabatic reactors are located 
under the hood the left the isothermal reactor, the back 
the unit. The instrument panel the middle contains most 


the pressure gauges, control valves and pressure regulators. 


wi 


Figure 9—Overall view the pilot plant. 


269 


INDUSTRIAL 


i 
) 


The stabilizer level control sight-glass and probe can seen 
the left this panel. the foreground are the two air purged 
control panels housing the electrical equipment. 


Material Balance 


The catalyst performance measured the yield and 
uality the desired products this case hydrogen and 
reformate. 

Usually there one feed stream, and three product streams 
make gas, debutanizer gas and stabilized reformate. The recycle 
gas does not appear the material balance, but its composition 
the same that the make gas. These streams are all 
measured, and the product streams are analyzed. 


The composition the make gas and debutanizer gas 
determined mass spectroscopy, and expressed mole 
volume The composition the stabilized reformate 
measured Vapor Partition Chromatography and 
expressed weight Other inspections the reformate 
are the specific gravity, ASTM distillation curve, vapor pressure 
and octane number. 

IBM 650 digital computer used calculate the yields 
the products, and the quality the C;+ reformate. The 
latter cannot obtained direct measurements, because there 
not perfect separation the stabilizer. 

Table summarizes the input data the computer. They 
consist the analyses the individual streams, their quantities, 
conditions which they were measured and the operating 
conditions. 

Table shows the items calculated the computer. The 
yield the reformate calculated adding the amount 
pentanes and heavier the make gas and debutanizer gas, 
and subtracting the amount butanes the reformate from the 
yield stabilized reformate. checked summing the 
individual component yields calculated from the analy ses. 
The RVP reformate calculated adding butanes 
the reformate. 

The last item this figure the number moles hydro- 
cracked may require some explanation. Because the large 
number components the naphtha, and the variety possible 
reactions, very difficult, not impossible, calculate the 
extent any one reaction. But every time molecule cracked, 
two molecules are produced. Therefore, the difference between 
the moles hydrocarbons the products and the feed repre- 
sents the net number individual cracking steps, C-C bonds 
broken. 

Table shows the basis this calculation. amount 
aromatics formation can expressed the number moles 
hydrogen produced adding the number moles hydro- 
cracking the measured production drogen. The extent 
naphthene dehydrogenation and dehydrocyclization 
cannot differentiated these data, but this can done from 
analysis feed and products. 

The results the individual calculations are printed out, 
and summary table also printed the same time, suitable 
for inclusion reports shown Tables and 5.* The IBM 
program has been real time-saver, cutting the computing time 
from about man-hours/weight balance hand calculation 
few minutes machine time, thus speeding the evaluation 
the experimental results. 


Unit Performance 


The pilot plant has been service for about years, and 
has performed quite satisfactorily. has shown good repro- 
ducibility, and the results have agreed well with actual plant 
performance. Normally, the overall material balances have been 


and this paper have been deposited Document No. 
6877 with the ADI Auxiliary Publications Project, Photoduplication 
Service, Library of Congress, Washington 25, D.C. A copy may be 
obtained by citing the Document No. and by remitting $1.25 for photo- 
prints, or $1.25 for 35 mm. microfilm. Advance payment is required. 
Make cheques or money orders payable to: Chief, Photoduplication 
Service, Library Congress. 
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TABLE 


Feed: 
Volume and blowcase temperature 
Gravity and distillation curve 


Gases: 
Wet test meter volume and temperature 
Sample bomb volume and pressure 
Analyses and atmospheric pressure 


Reformate: 
Weight reformate 
Gravity and distillation curve 
Octane number and vapor pressure 
Analysis 


Miscellaneous: 
Weight catalyst 
Hours weight balance 
Reactor pressure 
Recycle gas rate 


TABLE 
IBM 650 CALCULATED RESULTS 


Weight Balance: 
Molecular weight and wt. feed the three streams 
Material loss 


Yield Individual Components: 
Wt. Vol. and Mole feed 


Reformate: 
Yield Wt. Vol. and Mole feed 
Octane number, density and vapor pressure 


Yield Vol. feed 
Octane number 


Miscellaneous: 

Recycle gas rate MSCF/B 

Feed rate W/W. hr. 

Hydrogen partial pressure 

Moles hydrocracked/100 moles feed 


TABLE 


METHOD FOR CALCULATION TOTAL PRODUCED 


No. Moles Increase 


Aromatic Formation Hydrogen Hydrocarbon 

Hydrocracking 


Data required from material balance 


within 98-102%, and frequently within 99-101%. The main 
problem this respect has been the maintenance leak-proof 
joints with frequent heating and cooling. 

This pilot plant not restricted catalytic reforming, 
and could used study any reaction involving fixed bed 
catalyst and liquid gaseous feeds elevated temperatures and 
pressures. Minor modifications may required, such steam 
tracing, high melting point materials were used, such 
waxes. 
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CHEMICAL ENGINEERS ESTABLISH AWARDS 


Commencing this autumn, the Chemical Engineering Division the will 
present two first for the best paper published The Canadian Journal 
Chemical Engineering during the year immediately preceding the Fall Conference, 
and the second for the best oral presentation technical paper the Annual 
Conference (Chemical Engineering Session) and the Fall Conference the Chemical 
Engineering Division. The awards will take the form engraved certificate and 
silver tankard with the proper inscription. 


For the first award, the winning paper will selected from the October issue 
through August issue C.J.Ch.E. restrictions are placed the nature 
che paper (whether research, industrial development review type), nor the 
nationality, and affiliation the author authors. case joint authorship, each 
the authors would receive identical certificates and tankards. Three judges have 
been set represent academic circles, industry, and government, have fair 
representation the three major areas chemical engineering activity Canada. 


For the second award, oral presentation, the program chairman appoint 
three judges rate the various papers, only single sessions are held. These judges 
are represent industry, academic circles, and government, alternately research 
and development organizations which not properly fit the first two categories. 
uniform method scoring will used and the winner will announced the 
end the meeting. Formal presentation will then made the next Con- 
ference. simultaneous sessions are held will necessary appoint one more 
judges order that all sessions will screened carefully. 


q 
Ag 
q 
q : 
3 
: 
; 
: 
4 
: 
iF fi 
q 
J 
4 
, 
2 
a 
> 
Ps 
H 


q 
q 


